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ABSTRACT 


Carbon nanotube field ionization teehnology has the potential to make ion propulsion 
feasible for use in miero- and nano-satellites. To better understand the phenomenon and 
optimize the ion thruster design, it is useful to have an aeeurate model of the system. 
Numerieal modeling of large-seale eleetron bombardment ion engines is a relatively 
mature field, but modeling of field-ionized ion engines is in its infaney. A simpler eode 
may be appropriate for the early modeling stages of earbon nanotube field ionization 
teehnology; one sueh software paekage is spiffe. Spiffe is intended for modeling 
axisymmetrie radio frequeney guns, but it eontains all the eode neeessary for basie 
modeling of ion opties in a field-ionized ion thruster. 

This work analyzes the feasibility of spiffe software for use in modeling field- 
ionized ion thrusters. It also provides detailed proeedures for its use. In this work, 
spiffe is first verified to agree with theoretieal predietions of limits in a one¬ 
dimensional approximation using eleetrons. Two primary geometries and their boundary 
eonditions are investigated. The geometry is then varied to test the limits of the one¬ 
dimensional approximation. This was further altered to simulate singly-eharged Argon 
ions and verified against theoretieal one-dimensional limits. 

A supplemental user’s guide was developed to aid students with minima! 
programming experienee to quiekly beeome familiar with the methods used in spiffe 
and the impaet of program options on the results. A guide to quickly post-processing the 
data was also developed. 
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I. 


INTRODUCTION 


A. BACKGROUND 

Carbon nanotube field ionization technology has the potential to make ion 
propulsion feasible for use in micro- and nano-satellites. To better understand the 
technology and also to optimize the ion thruster design, it is useful to have an accurate 
model of the entire system. Numerical modeling of standard large-scale electron 
bombardment ion engines is a relatively mature field, but modeling of field-ionized ion 
engines is in its infancy. Current numerical modeling schemes are concerned with 
characteristics of a plasma, the ion source, at the plasma sheath interface. Utilization of 
an available code may be appropriate for the early modeling stages of carbon nanotube 
field ionization technology; one such software package isspiffe. Spiffeis intended 
for modeling axisymmetric RF guns, but it contains all the code necessary for basic 
modeling of the ion optics in a field-ionized thruster. 

B. OBJECTIVES 

Any numerical modeling scheme requires a thorough understanding of the 
applicability and the limits of the program. This logically begins with a study of the 
equations and boundary conditions utilized, and a comparison to known or theoretical 
results. This work analyzes the feasibility of spif f e software for use in modeling field- 
ionized ion thrusters. The research questions to be addressed have been: 

• How does spiffe differ from current steady-state numerical modeling 
schemes for ion optics? 

• Does spiffe accurately model a theoretical one-dimensional situation? 

• Can spif f e be altered to model positive ions? 

• What are the limitations of using spiffe to model such ions? 
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C. ORGANIZATION 

Chapter II provides a background on electrostatic thrusters, focusing on relevant 
equations and theoretical limits. Concepts specific to ion thruster design are introduced. 

Chapter III provides an overview of spiffe. This section is designed to give 
the user enough background on spiffe to understand the methodology presented in 
Chapter IV. A more detailed description of spiffe is provided as a supplemental 
user’s guide in Appendix A. Chapter III also answers two of the research questions. 
First, how does spiffe differ from current modeling schemes? Second, can spiffe be 
altered to model positive ions? 

Chapter IV outlines the research approach and provides a description of initial 
program set-up. Chapter V presents data from the initial set of simulations. A 
preliminary analysis of data from Chapter V was used to shape the follow-on simulations, 
the results of which are presented in Chapters VI and VII. 

The first part of the study focused on all commercially available modeling 
programs and their applicability in modeling carbon nanotube field-ionized thrusters. 
This is covered in Chapter II. This survey of available modeling programs led to 
selection of spiffe for further study. 

The second part of the study focused on familiarization of spiffe, program 
structure, method of coding, equations, and boundary conditions. The result of this study 
was the development of a supplemental user’s guide (Appendix A). This guide was 
developed to aid students with minimal programming experience in quickly becoming 
familiar with the methods used in spiffe and the impact of program options on the 
results. A guide to post-processing the data was also developed. Chapter III provides an 
abbreviated overview of spiffe and its methods, as well as the general set-up for the 
simulations. 
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Part three used a one-dimensional approximation to verify that spiffe agreed 
with standard results in literature. This part of the study explored various geometries and 
the impact of boundary conditions. This led to selection of a simplified geometry for 
further study. 

The simplified geometry was varied to test the applicability of the one¬ 
dimensional approximation. The code was also altered to simulate singly-charged Argon 
ions and again verified against theoretical one-dimensional limits. 
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II. BACKGROUND 


A. ELECTROSTATIC THRUSTERS 

Electrical propulsion systems are grouped into three major categories: 
electrothermal, electrostatic, and electromagnetic. Electrostatic thrusters rely on the 
Coulomb force to accelerate propellant ions. An ion thruster schematic is shown in 
Eigure 1. This model assumes one-dimensional ion flow in the ^-direction. The source 
electrode is located at x=0; the accelerating electrode is located at x=d. The potential 
difference between the source and accelerating electrode is V. The method of ion 
production is not addressed here; ions are assumed to exist at the source in a uniform 
distribution. The potential difference between the ion source and the accelerating 
electrode establishes an electric field. The electric field accelerates the ions from the 
source through the gaps in the accelerating electrode. The ions are neutralized as they 
exit the thruster to prevent charge build-up. 



Eigure 1. Schematic diagram of an ion thruster. Erom [1]. 


The one-dimensional electric field established by the electrodes is given by 



The force F that is experienced by an ion with charge q is 

F=qE (2) 

The particle of mass m will be accelerated according to Newton’s Second Eaw. 
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F 

a = — 
m 


(3) 


Assuming that the particles are emitted with zero velocity at the source, conservation of 
energy requires that the kinetic energy of the particle equals the electrical energy gained 
in the field. 


1 2 Tr 

—mv -qV 

2 


(4) 


The time required for the particle to traverse the interelectrode space (d) is given by 


t = 



(5) 


Combining the above equations, it can be shown that the total thrust generated by the 
accelerated particles is given by 


F = 



( 6 ) 


where 1 is the current of the ion beam. 

1. Space Charge 


Equation (6) suggests that the thrust can be increased without bound by simply 
increasing the current, but this is not possible. In addition to the external electric field 
established by the electrodes, each ion produces its own electric field. The electric field 
of the ions counteracts the external field created by the electrodes. This effect will be 
greatest where the ion density is greatest, at the source electrode. When the electric field 
from the ions counteracts the external field at the source, no further current may be drawn 
from the source. This is known as the space-charge-limited case and the resultant current 
density is known as the saturation current density. In the case of infinite parallel plates, 
the saturation current density (7) is given by the Child-Langmuir law. 


9 i m d^ 


(7) 


where is the permittivity of free space, q and m are the charge and mass of the particle, 

V is the accelerating voltage, and d is the distance between the source and the 
accelerating electrode. 
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For a given particle (fixed qlm), the saturation current density will depend only on 
the accelerating voltage (V) and the interelectrode spacing (d). 

Space-charge effects cause the electric field, potential, and ion velocity to vary 
from the ideal electrostatic case. Figure 2 illustrates the relationship between the ideal 
case and the space-charge-limited case. In the ideal case (no particles or negligible 
charge density), the electric field is constant and the potential increases linearly between 
the electrodes. In the space-charge-limited case, the electric field is nonhomogeneous 
and approaches zero at the source electrode. The potential distribution between the plates 
becomes nonlinear. 



k 




©■ 


Accderoting 

griC 
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Figure 2. One-dimensional ion acceleration for negligible charge density (—) and 

for space-charge limit (—). From [1]. 
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The Child-Langmuir law derives from the one-dimensional Poisson equation 


d^V _ p 
dx^ £q 


( 8 ) 


where prepresents the space charge density. A general solution is determined by solving 
the three-dimensional Poisson equation. 

d^V d^V d^V p 

Current density is related to space charge density by 

J = pv (10) 

where velocity (v) is given by Equation (4). 

The Poisson equation can be derived from Maxwell’s equations. The two 
relevant equations are Gauss’ law 


and Faraday’s law 


VD = p 


VxE = -— 
dt 


Assuming a homogeneous, isotropic, linear medium 


Gauss’ law becomes 


D = s^E 




( 11 ) 


( 12 ) 


(13) 


(14) 


In an electrostatic system, the magnetic field is zero. Faraday’s law becomes 

VxP = 0 (15) 

Equation (15) indicates that E can be written as the gradient of a scalar (potential). 

F = -VF (16) 

Combining equations (14) and (16) results in 

V-F = -VV = ^ (17) 


which is equivalent to equation (9). 
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2 . 


Beam Spreading 


The electric field from each ion also acts on surrounding ions. The interaction 
between ions is a repulsive force that causes beam spreading. The electric field (E) from 
a point charge is given by 


E = 


1 q . 

-yr 

Ans^ r 


(18) 


where So is the permittivity of free space and r is the radial distance from the point 
charge to the location of interest. 


3. Aspect Ratio 


For a circular ion beam of radius r, the current is given by 


I = Jnr^ 


= Jn 


f 

V 



Combining this relationship with Equations (6) and (7) gives 

^ 2;rs,V^ fDV 
9 yd) 


(19) 


( 20 ) 


The ratio D/d is called the aspect ratio (AR). It is a ratio of the effective beam diameter 
(D) to the interelectrode spacing (d) and uniquely relates the thrust to the square of the 
applied voltage [1]. Equation (20) suggests that the thrust may achieve large values 
simply by increasing the aspect ratio. In reality, as D increases, “potential gradients in 
the radial direction become so great that beam collimation by external electrodes is no 
longer feasible [2].” A practical limit on the aspect ratio for a single beam accelerated by 
an external electrode is on the order of unity. Modern ion thrusters have adopted a 
gridded accelerator structure to overcome this limitation. The ion beam that emerges 
from each aperture is termed a beamlet. The gridded structure allows each beamlet to 
approach the practical aspect ratio limit and maximize the thrust produced. The total 
thrust of the ion beam is a summation of the thrust produced by each beamlet. 
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4. 


Perveance 


Perveance is a measure of performance of the thruster and is based only on the 
electric parameters of the system. Perveance is defined as 


P = 



and is related to the aspect ratio by 


P = 


^^0 
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( 21 ) 


( 22 ) 


Perveance is essentially a measure of the current-producing capability of the 
system. It is useful because it is definable for any accelerator system, whereas aspect 
ratio is restricted to the plane diode case [1]. Given the perveance, an equivalent plane 
diode may be developed for more complex systems. 


5. Ion Optics 

All electrostatic engines consist of an ion source, a means of accelerating those 
ions, and a means of neutralizing the exhaust. Ion optics refers to the accelerating 
assembly. 

In designing the ion optics, it is essential to determine the aspect ratio that results 
in the optimum perveance for a given beam current and voltage. Alternatively, the design 
challenge may be to determine the optimum current to voltage ratio for a given aspect 
ratio. In either case, operating at other than the optimum perveance results in beam 
impingement on the grids, which both damages the grids and reduces the number of ions 
available to provide thrust. Figure 3 illustrates under-perveance (a), optimum perveance 
(b), and over-perveance (c). Over-perveance results in impingement on the accelerating 
grids; under-perveance results in crossover, overlap of neighboring beamlets, and 
impingement on downstream surfaces. Optimum perveance results in ideal focusing of 
the beam without impingement or crossover. 
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Figures. Comparison of perveance conditions. From [3]. 

B. ION SOURCES 

Ions may be produced by several methods: photo-ionization, electron 
bombardment, field ionization, surface ionization, and thermionic emission. Electron 
bombardment is a well-established technology [3] [4] and several numerical models have 
been developed to aid in the design of such thrusters. Field ionization thrusters are under 
development. Numerical models associated with field ionization typically deal with non¬ 
propulsion applications, such as mass spectrometry. 

1. Electron Bombardment 

In an electron bombardment thruster, ionization is accomplished as electrons 





11 







































































collide with gaseous propellant atoms in the ionization ehamber. Ionization of the 
propellant gas ereates both ions and eleetrons, forming a plasma. A simplified sehematie 
of an eleetron bombardment thruster is shown in Figure 4. Eleetrons are attraeted to the 
anode but are foreed by the magnetie field to spiral in the chamber, eausing numerous 
eollisions and more ionizations. If an eleetron internets with an ion, another ionization 
may oeeur, thereby produeing a doubly-eharged ion. Doubly-eharged ions have a 
different charge-to-mass ratio than the singly-eharged ions and thus a different trajectory 
through the ion opties and are eonsidered as losses. 



Figure 4. Simplified sehematie of an eleetron bombardment thruster. From [4]. 

Eleetrons are removed from the ehamber by the radial electric field. The axial 
eleetric field moves the ions toward the accelerator grids. Ions must first pass through the 
screen grid before being aceelerated by the accelerating grid. The purpose of the sereen 
grid is to protect the accelerating grid from off-axis ions and provide some confinement 
of the eleetrons in the plasma [3]. 

The presenee of plasma in these engines has a signifieant impact on the design of 

the ion opties. Figure 5 illustrates the various ways in which the shape of the plasma 

sheath ean serve to shape the ion beam. Each case represents a ehange in the accelerating 
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voltage. When there is no accelerating voltage applied (a), the plasma extends into the 
accelerating region and forms a sheath at each electrode. As the voltage is increased, the 
plasma sheath is pushed back toward the discharge chamber (b) and the plasma boundary 
takes on a convex shape (c). As the voltage is increased, the plasma boundary recedes 
further into the discharge chamber and takes on a concave shape (d). The concave shape 
is desirable because it provides initial beam shaping. 


(o) NO ACCEL VOLTAGE 




(b) LOW ACCEL VOLTAGE 




.ION 

SHEATH 


PLASMA 

BOUNDARY 


(dl NORMAL OPERATING ACCEL VOLTAGE 


Figures. Effect of plasma sheath on initial ion trajectories. From [5]. 

While the shape of the plasma sheath is not directly at the disposal of the 
designer, the ion trajectories and source shape are determined by the accelerating voltage 
and current density. The variation in the shape of the plasma sheath due to changing 
operating conditions must be considered in the electrode design. Electron bombardment 
thrusters utilize a tapered screen aperture so that the ion current and trajectories are as 
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insensitive as possible to the shape of the plasma sheath [6]. Figure 6 illustrates 
variations in the plasma sheath and the benefit of a tapered sereen eleetrode for different 
operating eonditions. 
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Figure 6. Illustration of the effeet of a tapered sereen eleetrode. From [5]. 

2. Field Ionization 

Field ionization occurs when an electron is stripped from a particle by a strong 
electric field. It is essentially field emission in reverse, and many parallels can be drawn 
between field ionization and surface ionization [6]. A schematic of a carbon nanotube 
(CNT) field-ionized thruster is shown in Figure 7. Propellant flows through a porous 
surface and past the CNTs. The external electric field is amplified at the tips of the 
CNTs, becoming strong enough to cause the valence electron to tunnel from a propellant 
atom into a CNT. The electron is absorbed by the metallic surface, leaving only the 
propellant ion to flow through the engine to provide thrust. This is different from a 
plasma source, from which electrons and neutral particles may escape to impinge on the 


14 
































































ion optics. A planar arrangement of CNTs makes it appropriate to model the ion souree 
as planar, with no concern for the shape of a plasma sheath. 


Propellant 


CNTs Accelerator Grids 



Figure 7. Schematic of a CNT field-ionized thruster (not to scale). From [7]. 

Field ionization is superior to electron bombardment in that it produces only 
singly-charged ions [6]. This allows the ion optics to be optimized for a given charge to 
mass ratio. In an electron bombardment thruster, a doubly-charged ion will follow a 
different trajectory and may impinge on the ion optics. Numerical simulations must 
account for the doubly-charged ions when determining particle trajectories. This is not a 
concern with field ionization sources. The field ionization chamber designer need only 
be concerned with maximizing the propellant mass utilization factor, which measures the 
number of ions that can be delivered to the accelerator [8]. 
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III. INTRODUCTION TO SPIFFE 


A. OVERVIEW 

Spiffe (Space charge and /ntegration Forces For Flectrons) is a 2.5-D fully 
electromagnetic particle-in-cell (PIC) code for radiofrequency electron gun (RF gun) 
simulation and similar systems with cylindrical symmetry. The program is inherently 
unsteady and uses Maxwell’s equations and the Lorentz force to compute the time- 
evolution of fields and particles in a metal cavity. It has been used at NPS to study 
transmission of intense electron beams through apertures [9] and suppression of current 
fluctuations in an intense electron beam [10]. Spiffeis open source software written in 
C programming language. 

B. PARTICE-IN-CELL (PIC) METHOD 

PIC codes assume that one simulation particle consists of many physical particles. 
The simulation particle is termed a macroparticle or superparticle. It exists only as a 
point charge; it has no volume. The computer code solves equations of motion and 
calculates macroquantities (e.g., density, current density) from the position and velocity 
of the macroparticles. The charge-to-mass ratio of a macroparticle is the same as the 
charge to mass ratio of the physical particle, which means that the macroparticle follows 
the same trajectory as the physical particle. The generic scheme of a PIC simulation is 
illustrated in Figure 8. 


Input/Outpirt 



Figure 8. Generic order of operations for a PIC scheme. From [II]. 
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1 . 


Particle Mover 


The portion of the code that solves the equations of motion and follows the 
trajectory of each macroparticle is termed the particle mover. One method of coding the 
particle mover, and the method employed in spif f e, is the leap-frog method. The leap¬ 
frog method is an explicit solver, meaning that it depends on forces from the previous 
time step. Explicit solvers are simpler and faster than implicit solvers, but their stability 
requires a smaller time step [11]. The simulation is divided into discrete time moments 
(At) and the physical quantities (A) are calculated according to: 

t —^ tu kAt 

(23) 

A(t)^A,=A(t = t,) 

where k=0,l,2..., t is the time, to the initial moment, and A denotes any physical quantity. 
The leap-frog method calculates some of the physical quantities (B) between time steps. 

h+\n =to + {k + l/ 2) At 

(24) 

B(t ) ^ = h+V2 ) 

For example, the particle position may be calculated at each time step and the 
particle velocity between time steps. In this example, particle position is calculated using 
the velocity from the previous half-step. The particle position is updated and the new 
particle velocity is calculated from the new position. 

2. Particle and Field Weighting 

Particle weighting and field weighting are terms used to describe the calculation 
of macroquantities and field data inside the simulation region. The simulation space is 
divided into a grid and the desired data is calculated at grid points or interpolated 
between grid points. 

C. SPIFFE METHODS 

A robust description of the equations and methods used in spif fe is provided in 
[12]. This section is intended as a brief overview of the basic methods used in spif f e. 
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1 . 


Particle and Field Propagation 


This section summarizes the methods described in [12]. In spiffe, the particle 
mover uses position (z, r) and momentum (p^, pr). The position data is calculated at each 
time step (t=to+nAt) and advanced to the next time step (t=to+(n+l )At) according to 


z > z \ At 


r > r I At 


r 

Ei£ 

r 


(25) 


where c is the speed of light, /is the Lorentz factor, and p=Py. Momentum is given in 
terms of electron mass and the speed of light {niec). 

The electric field is calculated at each time step; the magnetic field is calculated 
between time steps. At each time step, the electric fields are calculated using the 
magnetic field and current values from the previous half-step. The magnetic fields are 
then advanced to the next half-step using the just-calculated electric fields. 


Particle-induced currents are needed to calculate the space-charge effect on the 
electric field, and these currents are computed from particle velocities. Particle position 
is stepped with the electric field and particle momentum is stepped with the magnetic 
field (between time steps). 


Pz^Pz+^^ 
P,^ Pr + At 
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2. Space Charge 

Particle-induced fields (i.e., space charge) are calculated using equivalent current 
density. Equivalent current density is calculated from the charges, momenta, and 
positions of the macroparticles. Spiffe locates the four grid points that surround the 
particle and allocates the charge between the upper and lower pairs using linear 

weighting. The current density is calculated from the charge and velocity of the 
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macroparticle and the volume of the cell surrounding the grid point. The resultant current 
density is split between the left and right pairs using linear weighting [12]. 

Spiffe allows the user to run simulations with or without space charge 
included. With the space charge off, the simulation provides results that are consistent 
with basic electromagnetics. 

3. Poisson Correction 

The field propagation routine uses two of Maxwell’s equations: Ampere’s law 
and Faraday’s law. It does not use the divergence equation. Accumulated errors in the 
field propagation routine lead to non-zero divergence. Spiffe corrects for this by 
solving the Poisson equation at user-defined intervals during the simulation and clearing 
any fictitious charge build-up. 

4. Lost Particles 

Spiffe is strictly an axi-symmetric code. Particles that impact a metal surface or 
outer boundary are lost. Particles that are lost to the centerline (r=0) are reflected back 
into the simulation region. 

5. Boundary Conditions 

Reference [13] indicates that specification of the potential on a closed surface 
defines a Dirichlet boundary condition, while specification of the electric field on the 
surface defines a Neumann boundary condition. In spiffe, the user is only able to 
specify the potential of boundaries, which implies that Dirichlet is the only possible 
boundary condition. However, [12] uses Maxwell’s equations and properties of even and 
odd functions to show how both Neumann and Dirichlet boundary conditions are 
implemented in spiffe. When a Neumann boundary condition is specified in 
spiffe, the electric field is perpendicular to the boundary. Conversely, the Dirichlet 
boundary condition requires that the electric field be parallel to the boundary. 
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6 . 


Geometry Definition 


The user defines the geometry using a series of points, which are connected by 
lines or arcs. Surfaces inside the cavity are assumed to be metal; there is no provision for 
changing the type of material. The source electrode may be placed anywhere; it does not 
need to correspond to physical boundary. 

7. Modifications for Ions 

Spiffe provides an option for increasing the mass of the particle to simulate ions. 
This is accomplished through use of the “stiffness” option. There is no provision for 
changing the charge of each particle. This implies that only singly-charged ions can be 
simulated. 
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IV. METHODOLOGY 


A. RESEARCH APPROACH 

1. Analysis of Simulation Tools 

Spi f f e was installed on a Windows-based laptop computer. The first part of the 
study involved a brief investigation of various electrode configurations to select which of 
them was most suitable for the problem at hand. 

2. One-Dimensional Approximation 

In order to test the simulation performance, a simple one-dimensional scenario 
was created to verify basic electrostatics (space charge off) and agreement with the 
Child-Langmuir law (space charge on). With the space charge off, the delivered current 
density should match the requested current density. This verifies that the code is working 
properly and that post-processing of the data was done correctly. With the space charge 
on, the delivered current density should be limited to the saturation current density given 
by the Child-Langmuir law. 

In order to approximate a one-dimensional situation, a parallel plate geometry 
with an aspect ratio of 100 was created. The cathode radius was based on anticipated 
modeling of the setup described in [8]. 

The theoretical saturation current density is calculated from the perveance, which 
is determined using the aspect ratio. First, the perveance is calculated using Equation 
(22). Then the current is calculated using Equation (21). The current is converted to 
current density using the cathode radius. 

As spiffe was designed for electrons, all simulations were done with electrons 
as the particle being accelerated. Chapter IX investigates the possibility of using 
sp i f f e to simulate ions. 
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3. Geometry Variations 


After the model was verified to agree with theoretical predictions, the geometry 
was varied to test the extent of the applicability of theoretical predictions. First, the 
aspect ratio was varied by increasing and decreasing the interelectrode spacing. Next, the 
aspect ratio was maintained at 100 while the geometry was scaled up. 

4. Ion Variations 

In order to investigate the ability of spiffe to model ion flow, the simple box 
geometry was again utilized. The initial set of simulations (AR=100, varying current 
density) was repeated for both space charge on and off. 

B. INITIAL CONDITIONS 

A spiffe run has several variables for the user to select, including grid size, 
cathode placement, initial momentum, simulation time. These can be overwhelming 
when trying to select the “right” initial conditions for a simulation. The following initial 
conditions were selected arbitrarily. 

1. Cathode Placement 

In all simulations, the physical cathode was drawn at z=0. Particles were released 
slightly to the right of the cathode. 

2. Initial Momentum 

Spiffe allows the user to specify an initial momentum of the particles when 
they are released from the cathode. Initial momentum was set to zero for all simulations. 

3. Grid Size 

Proper grid sizing is essential to balancing accuracy and reasonable simulation 
time. A mesh that is too coarse will provide quick but inaccurate results. A mesh that is 
too fine will be prohibitive in terms of CPU memory and simulation time. The grid lines 
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were selected to provide a resolution approximately one to two orders of magnitude less 
than the smallest dimension of geometry. This resolution provided the best balance of 
accuracy and reasonable run time. Grid sizing is addressed in detail in Appendix A. 

4. Time Step 

Time step is directly related to grid size. If the number of grid lines doubles, the 
time step required for stability is halved. As the size of the time step decreases, spiff e 
must perform more calculations. This results in a longer run time, which can exceed 
twelve hours in some cases. The grid size, simulation time, and time step must be chosen 
carefully to balance accuracy with run time. 

Spiffe provides feedback on the minimum time step required for stability. In 
all simulations, a time step was selected to meet the stability requirement. 

5. Current Density 

Spiffe allows the user to specify a constant current density or to load a current 
density profile (e.g., half-gaussian ramp-up). A constant current density was used for all 
simulations. Current density was the primary variable in the simulations. 

This report frequently compares current density requested versus current density 
delivered. The requested current density is the value that is specified by the user in the 
spiffe code. The current density delivered is the output from the simulation. It is 
calculated by recording the total charge to pass through a screen over a given time period. 

6. Macroparticles 

Spiffe allows the user to specify macroparticles per step or electrons per 
macroparticle. These options are disussed further in Appendix A. In all simulations, 
eight macroparticles per step was used. 
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7. 


Simulation Time 


Simulation time refers to the time specified in the spiffe code; run time refers 
to the actual time that the simulation took to finish. Selection of the simulation time is 
critical for ensuring enough time has passed to allow the simulation to reach steady-state, 
but not so long as to lead to excessive run times. 

Equation (5) provides the flight time for a particle. The simulation time was 
selected by first determining the flight time of a particle. The flight time was increased 
by a factor of 15 and rounded up slightly to give the simulation time. This method 
provided a reasonable balance of accuracy and run time. 

8. Boundary Conditions 

The default boundary conditions were used for all boundaries. These apply only 
to the cavity boundaries and not to any of the surfaces inside the cavity. The upper, left, 
and right boundaries are Neumann; the lower is Dirichlet. 

• Upper: Neumann, Ej, = 0 at r = r^ax 

• Left: Neumann, Er = 0 at z = Zmin 

• Right: Neumann, Er = 0 at z = Zmax 

• Lower (centerline): Dirichlet, Er = 0 atr = 0 

9. Accelerating Voltage 

Early simulations used an arbitrary voltage that resulted in electric field strength 
of 1 MV/m. As geometries were varied, the voltage was adjusted to maintain the electric 
field strength at 1 MV/m. 

C. DATA CAPTURE 

The primary method of extracting data from spiffe was the use of diagnostic 

screens. Screens capture location, charge, momentum, and position data for every 
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macroparticle in the beam. Three screens were placed in the beam path to capture 
particle data, one near the cathode, one near the anode, and one in the middle of the gap. 
The screens extend from r=0 to r=rmax- The current density delivered was calculated 
from the data obtained from the screens. Plots included in this report use the current 
density from the middle screen. 

D. BASELINE GEOMETRIES 


I. Simple Box 


The simplest geometry is a box that is 0.1 mm wide (z) and 10 mm tall (r). The 
cathode is located at 0 < r < 5 mm. 



2(5 mm) 
0.1 mm 


= 100 


(27) 


When the simulation starts, electrons are accelerated from the cathode (0 V) to the 


anode (100 V). The cathode and anode are simply the left and right borders of the box, as 
shown in Figure 9. 





Figure 9. Electron beam acceleration in the simple box geometry. 


a. Boundary Conditions 

The left and right boundaries may be set to either Neumann or Dirichlet, 
but the upper is limited to Neumann and the lower (centerline) is limited to Dirchlet. In 
the simple box geometry, the left and right boundary conditions must be set to Neumann. 
If they were set to Dirichlet, then would be zero at the left and right boundaries, and 
there would be no field to accelerate the electrons across the gap. 
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b. Screens 

Three screens are placed in the beam path to capture particle data. The 
screens are located at 0.01 mm, 0.05 mm, and 0.09 mm. 

2. Parallel Plate 

The parallel plate geometry uses the same interelectrode spacing and cathode 
radius as the simple box. The plates are 0.01 mm thick and are separated by 0.1 mm. 
The cathode is assigned a potential of 0 V; the anode 100 V. The geometry and beam are 
shown in Figure 10. 



z (m) z (m) z (m) 

Figure 10. Electron beam acceleration in the parallel plate geometry. 


a. Boundary Conditions 

The boundary conditions apply to the borders of the cavity. As with the 
box geometry, they are set to their default values. Any boundaries drawn inside the 
simulation region are assumed to be metal surfaces. As such, the transverse electric field 
at the surface of the plates is zero. 

b. Screens 

The screens are placed at the same location as in the box geometry: 0.01 
mm, 0.05 mm, and 0.09 mm. 
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V. ONE-DIMENSIONAL APPROXIMATION 


A. BOX - SPACE CHARGE OFF 
1. Current Density 

With the space charge off, spiffe returns the requested current density (Figure 
11). For example, if the code specified 200 kA/m , spiffe delivered approximately 
200 kA/m . This verifies that the code is working properly and that post-processing of 
the data was done correctly. 



Current Density Requested (A/m^) 

Figure 11. Current density requested vs. delivered, box geometry, space charge off. 

2. Electric Field 

With the space charge off, the electric field between the cathode and anode should 
be constant. 

E = = 1x10® y/m (28) 

^ 0.1 mm 


The electric field at the centerline {r=0) as a function of z is shown in Figure 12. 
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Each plot represents a time snapshot at intervals of 0.05 ns and all are plotted to the same 
scale. With the space charge off, the magnitude of the electric field remains constant at 1 
MV/m. 



The plot of the equipotential lines, Figure 13, indicates a constant field for all z 

and r. 
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Figure 13. Equipotential plot, Ez as a function of z and r, Jreq=300 kA/m^, space 

charge off. 


3. Screen Data 

Data from the screens is contained in Appendix B. This data shows that the 
average longitudinal and radial momenturns do not vary with requested current density. 
The maximum radial coordinate of the particles does not exceed 5 mm and is constant for 
all values of requested current density. 

B. BOX - SPACE CHARGE ON 


1. Current Density 


With the space charge on, the simulation should be limited by saturation current 
density. The expected saturation current density is calculated using Equations (12) and 
(13). 


P = ^J%ARf = 


;r(8.85xlO-'') 12(1.602x10-“') 
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-31 
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(31) 




Figure 14 shows that current density for the described configuration, simulated in 

spiffe, levels off at a slightly greater than expected saturation current density, 

2 

approximately 270 kA/m . 



Current Density Requested (kA/m^ 


Figure 14. Current density requested vs. delivered, box geometry, space charge on. 

2. Electric Field 

Electric field snapshots are shown in Figure 15. Figure 16 is a magnified plot of 
the seventh plot. Each plot represents a time snapshot at intervals of 0.05 ns; the first plot 
shows the electric field before particles are released. The requested current density is 300 
kA/m . The expectation is that the electric field will be zero at the cathode, as illustrated 
in the generic case shown in Figure 2. The electric field does not go to zero at the 
cathode, but it decreases significantly and the field takes on the general shape of the 
expected distribution. 
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Figure 15. Ez at r=0 for Jreq=300 kA/m^, space charge on. Plots are taken at time 
intervals of 0.05 ns. Figure 16 is provided for scale. 
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Figure 16. Ez at r=0 for Jreq=300 kA/m^, space charge on (snapshot 7 from Figure 15). 

Electric field contour plots are shown in Figure 17. These show significant 
interaction between the electric field of the particles and that of the electrodes. 
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Figure 17. Equipotential plots, Ez as a function of z and r, Jreq=300 kA/m^, space 
charge on. Plots are taken at time intervals of 0.05 ns. 

3. Screen Data 

Data from the screens is contained in Appendix B. With the space charge on, the 
average longitudinal and radial momentums vary with requested current density. The 
maximum radial coordinate of the particles exceeds 5 mm, indicating beam spreading. 

C. PARALLEL PLATE - SPACE CHARGE OFF 
1. Current Density 

With the space charge off, spiffe returns the requested current density (Eigure 
18). Eor example, if the code specified 200 kA/m , spiffe delivered approximately 
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2 

200 kA/m . This verifies that the eode is working properly and that post-proeessing of 
the data was done eorrectly. This also provides a eomparison to the box geometry. 



Current Density Requested (kA/m^ 


Figure 18. 

2 . 


Current density requested vs. delivered, parallel plate, space charge off. 

Electric Field 


The electric field as a function of z at r=0 is shown in Figure 21. The electric 
field at the centerline is constant (£’z=l MV/m) for all points between the plates, with the 
exception of the two points closest to the electrodes (£'z=0). These points did not appear 
in the box geometry. These points are located on the grid, which does not correspond 
to the grid that defines the physical boundaries (see Appendix A). 
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Figure 19. as a function of z at r=0. Plots are taken at time intervals of 0.05 ns. 

Figure 20 is provided for scale. 



z (m) 


Figure 20. as a function of z at r=0 (snapshot 7 from Figure 19). 


























The electric field contour plot is shown in Figure 21. (This plot only shows the 
contours from z=0 mm to z=0.1 mm and r=0 mm to r=10 mm in order to match the 
contour plot from the box geometry.) This plot indicates a difference from the box 
geometry. Each plot uses 20 contour lines, but the parallel plates do not show the 
regularly spaced equipotential lines seen in Figure 13. 
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Figure 21. Ez as a function of z and r, parallel plate, space charge off. 

3. Screen Data 

Data from the screens is contained in Appendix B. Despite the differences in the 
electric field, the screens return the same data as the box geometry. 

D. PARALLEL PLATE - SPACE CHARGE ON 

1. Current Density 

The expected saturation current density is the same as that calculated in Equation 
(23). Figure 22 shows that spiffe levels off at a saturation current density less than 
expected, approximately 190 kA/m^. 
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Figure 22. Current density requested vs. delivered, parallel plate, space charge on. 

2, Electric Field 


The electric field at the centerline (r=0) as a function of z is shown in Figure 23. 
Each plot represents a time snapshot at intervals of 0.05 ns; the first plot shows the 
electric field before particles are released. The requested current density is 300 kA/m . 
The electric field again takes on the distribution from Figure 2, with the exception of the 
two points near the electrodes. These points remain at zero, as they were when the space 
charge was off (Figure 19). Figure 24 shows a zoomed plot of the seventh plot of Figure 
23. 
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Figure 23. as a function of z at r=0, parallel plate, space charge on. Plots are taken 

at time intervals of 0.05 ns. Figure 24 is provided for scale. 
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Figure 24. Ez as a function of z at r=0, parallel plate, space charge on (snapshot 7 
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from Figure 23). 
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The equipotential plots (Figure 25) show differences from the box geometry. In 
general, the equipotential lines from the box geometry are more evenly spaced. 



I [mj I (m) 


Figure 25. Ez as a function of z and r, parallel plate geometry, space charge on. Plots 

are taken at time intervals of 0.05 ns. 


3. Screen Data 

Data from the screens is contained in Appendix B. This data provides insight 
regarding the differences between the box and parallel plate geometries. 
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VI. PARALLEL PLATE VARIATIONS 


A. OVERVIEW 

The parallel plate and box geometries were expected to provide similar, if not the 
same, results. Since this was not the case, several additional geometry variations were 

'y 

investigated. Each variation used a requested current density of 300 kA/m . This value 
was selected because it is above the saturation current density. Values below the 
saturation current density have minimal impact on the results. 

B. REMOVAL OF ONE PLATE 

Removal of either the left or right plate results in a saturation current density that 
falls between the saturation current densities calculated in the case of the box geometry 
and the full parallel plate geometry. 

1. Removal of Left Plate 

When the left plate is removed, the left boundary becomes the cathode. The 
boundary condition at z=0 is Neumann (Er=0). The gap is still 0.1 mm. The potential of 
the cavity boundary is set to 0 V with the default boundary conditions. The remaining 
plate is set to 100 V and the potential of the cavity boundaries is set to 0 V. 

0.01 0 

0.008 
^ 0.006 
0.004 

O.n 

0.||''I 

z (m) z (m) z (m) 

Figure 26. Left plate removed, electron beam acceleration. 
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With the space charge off, the results (Table 1) are consistent with the parallel 
plate and box geometries. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (tOgC) 

Pr,avg (tOeC) 

rmax(mm) 

300 

0 

299.4907 

6.256E-03 

0 

4.999995 


1 

299.4482 

1.399E-02 

0 

4.999995 


2 

299.5006 

1.877E-02 

0 

4.999995 


Table 1. 

Left plate removed, space charge off. 



With the space charge on, the current delivered falls between the box geometry 
and full parallel plate geometry. The results are listed in Table 2. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (tneC) 

Pr,avg (tneC) 

rmax(mm) 

300 

0 

218.2687 

3.239E-03 

3.299E-06 

5.001634 


1 

215.8934 

1.132E-02 

1.688E-05 

5.011986 


2 

215.3945 

1.717E-02 

1.626E-05 

5.021632 


Table 2. Left plate removed, space charge on. 


2. Removal of Right Plate 


When the right plate is removed, the cathode is the left plate and the right 
boundary becomes the anode (Figure 27). The boundary condition at z=Zmax=0.1 mm is 
Neumann (Er=0). The potential at the left, top, and bottom boundaries is 0 V; the right 
boundary is set to 100 V. 




z (m) 

Figure 27. Right plate removed, electron beam acceleration. 
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With the space charge off, the results (Table 3) are consistent with the previous 
geometries. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (tOeC) 

Pr,avg (tOeC) 

rmax(mm) 

300 

0 

299.4907 

6.256E-03 

0 

4.999995 


1 

299.4482 

1.399E-02 

0 

4.999995 


2 

299.5006 

1.877E-02 

0 

4.999995 


Table 3. Right plate removed, space charge off. 


With the space charge on, the current density delivered falls between the box 
geometry and full parallel plate geometry, and is the same (to 4 significant figures) as the 
left-plate- removed scenario. The results are listed in Table 4. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (tneC) 

Pr,avg (tneC) 

rmax(mm) 

300 

0 

218.2281 

3.239E-03 

3.611E-06 

5.001601 


1 

215.9259 

1.132E-02 

1.716E-05 

5.011982 


2 

215.3820 

1.717E-02 

1.648E-05 

5.021624 


Table 4. Right plate removed, space charge on. 


C. PLATE THICKNESS VARIATIONS 

Returning to the original parallel plate geometry, the plate thickness was varied 
and the results compared to the baseline from Chapter V. Four scenarios were evaluated, 
two with increased thickness and two with decreased thickness. 

I. Increased Thickness 

The plate thickness was increased from 0.01 mm to 0.015 mm. The cavity was 
left at its original width, giving a 0.005 mm margin between the plates and boundary. 
This simulation failed to produce current. 

The simulation cavity was then widened by 0.005 mm on each side to give a 0.01 
mm margin between the plates and boundary. The grid resolution remained at 0.005 mm, 
providing a grid line on each boundary of the plate and one internal. 
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With the space charge off, the plate width had no effect. The electric field was 
consistent with the plots shown in the parallel plate section and the current density 
delivered was identical to the original geometry. The results are listed in Table 5. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz.avg (tOeC) 

Pr,avg (tOeC) 

rmax{mm) 

300 

0 

299.4907 

6.256E-03 

0 

4.999995 


1 

299.4482 

1.399E-02 

0 

4.999995 


2 

299.5006 

1.877E-02 

0 

4.999995 


Table 5. Plate thickness 0.015 mm, space charge off. 


With the space charge on, increasing the plate thickness had minimal effect on the 
current density. The results from the original geometry are shown in Table 6. The results 
from the thicker plates are shown in Table 7. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (tneC) 

Pr,avg (tneC) 

rmax(mm) 

300 

0 

187.3599 

3.112E-03 

4.479E-06 

5.002438 


1 

184.9992 

1.083E-02 

1.831E-05 

5.015919 


2 

185.8662 

1.641E-02 

1.830E-05 

5.027015 

Table 6. 

Plate thickness 0.01 mm (original geometry), space charge on 

Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (tneC) 

Pr,avg (tneC) 

rmax(mm) 

300 

0 

187.1261 

3.115E-03 

5.137E-06 

5.002460 


1 

185.0853 

1.083E-02 

1.839E-05 

5.015656 


2 

185.8173 

1.642E-02 

1.829E-05 

5.026623 


Table 7. Plate thickness 0.015 mm, space charge on. 


Taking it a step further, the plates were widened to 0.02 mm. This time the plates 
had two internal grid lines. The plate width had negligible impact on the current density, 
as shown in Table 8. 
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Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (rOeC) 

Pr,avg (rOeC) 

rmax(mm) 

300 

0 

187.0780 

3.115E-03 

5.737E-06 

5.002595 


1 

185.3433 

1.083E-02 

1.885E-05 

5.016203 


2 

185.9997 

1.641E-02 

1.857E-05 

5.027237 


Table 8. Plate thickness 0.02 mm, space charge on. 


2. Decreased Thickness 

The plate thickness was decreased from 0.01 mm to 0.005 mm. The grid 
resolution was left at 0.005 mm, providing a grid line on each border of the plate but 
none internal. With the space charge off, narrowing the plate had no effect, as shown by 
the data in Table 9. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (rneC) 

Pr,avg (rneC) 

rmax{mm) 

300 

0 

299.4907 

6.256E-03 

0 

4.999995 


1 

299.4482 

1.399E-02 

0 

4.999995 


2 

299.5006 

1.877E-02 

0 

4.999995 


Table 9. Plate thickness 0.005 mm, space charge off. 


With the space charge on, decreasing the plate thickness resulted in a decrease in 
current density of approximately 12%, as shown in Table 10. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (rneC) 

Pr,avg (rneC) 

rmax(mm) 

300 

0 

164.6370 

3.039E-03 

6.37E-06 

5.003461 


1 

164.6696 

1.051E-02 

1.96E-05 

5.018811 


2 

165.1050 

1.592E-02 

2.01E-05 

5.030431 


Table 10. Plate thickness 0.005 mm, space charge on. 


The mesh was adjusted to provide a grid resolution of 0.0025 mm longitudinally 
(keeping 0.005 mm radially). This provided one internal grid line for each plate. This 
resulted in a slight decrease in current (Table 11), which can be attributed to the change 
in grid resolution. 
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Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (fTleC) 

Pr,avg (fTleC) 

rmax(mm) 

300 

0 

159.4560 

3.363E-03 

4.37E-06 

5.003652 


1 

159.4374 

1.071E-03 

1.51E-05 

5.013048 


2 

159.8853 

1.602E-02 

1.50E-05 

5.022151 


Table 11. Plate thickness 0.005 mm, space charge on, 0.0025 mm longitudinal grid 

resolution. 

The plates were narrowed to 0.001 mm, resulting in another decrease in the 
current density. This data is shown in Table 12. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (fTleC) 

Pr,avg (fTleC) 

rmax(mm) 

300 

0 

144.9218 

3.563E-03 

2.943E-06 

5.001991 


1 

145.2404 

1.066E-02 

1.210E-05 

5.011490 


2 

145.6401 

1.581E-02 

1.230E-05 

5.019882 


Table 12. Plate thickness 0.001 mm, space charge on, 0.001 mm longitudinal grid 

resolution. 

There seems to be a critical point in plate width above which there is negligible 
impact on the current density. Tables 9-12 are summarized in Figure 28. 



Plate Thickness (mm) 


Figure 28. Variation of current density with plate thickness. 
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3. 


Cavity Size 


The margin between the plate and cavity border in the original parallel plate 
geometry is 0.01 mm. When the plate thickness was increased, the cavity was widened 
by the same amount to keep the margin between plate and cavity border at 0.01 mm. 
When the plate was widened without adjusting the cavity width, no current density was 
delivered. This led to an investigation of the impact of cavity size. The plate thickness 
was set at the original value (0.01 mm) and the cavity was widened to provide a margin 
of 0.02 mm between each plate and the cavity boundary. There was little impact on the 
current provided, as shown in Table 13. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (rOeC) 

Pr,avg (rOeC) 

rmax(mm) 

300 

0 

187.4077 

3.112E-03 

4.577E-06 

5.002438 


1 

184.9943 

1.083E-02 

1.825E-05 

5.015945 


2 

185.8647 

1.641E-02 

1.840E-05 

5.027093 


Table 13. Wide cavity (0.02 mm margin), original plate thickness, space charge on. 

D. PLATE RADIUS VARIATIONS 

This section investigates the impact of plate radius on the saturation current 
density. The beam radius remained at 5 mm while the plate radius was varied. 

1. Decreased Plate Radius 

All dimensions were held constant while the plate radius was decreased from 10 
mm to 7.5 mm. The beam radius remained at 5 mm. The reduction in plate radius 
caused a reduction in current density. The results from the original geometry are shown 
in Table 14. Table 15 lists the results from the decreased plate radius. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (f^eC) 

Pr,avg (f^eC) 

rmax(mm) 

300 

0 

187.3599 

3.112E-03 

4.479E-06 

5.002438 


1 

184.9992 

1.083E-02 

1.831E-05 

5.015919 


2 

185.8662 

1.641E-02 

1.830E-05 

5.027015 


Table 14. Original radius (10 mm), space charge on. 
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Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (f^eC) 

Pr,avg (f^eC) 

rmax(mm) 

300 

0 

147.7252 

2.970E-03 

4.777E-06 

5.008343 


1 

148.2254 

1.024E-02 

1.909E-05 

5.026992 


2 

148.1635 

1.549E-02 

1.993E-05 

5.040493 


Table 15. Decreased radius (7.5 mm), space charge on. 


The radius/height of the cavity was also varied. As with the thickness variation, 
this showed negligible impact. The results, listed in Table 16, are similar to the results 
listed in Table 15. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (rnqC) 

Pr,avg (rnqC) 

rmax(mm) 

300 

0 

148.7674 

2.985E-03 

3.294E-06 

5.006763 


1 

146.9948 

1.028E-02 

1.779E-05 

5.023468 


2 

147.4970 

1.551E-02 

1.919E-05 

5.033166 


Table 16. Decreased radius (7.5 mm), decreased cavity height (8.5 mm). 


2. Increased Radius 

The radius of the plates was increased to 20 mm; the results are listed in Table 17. 
The current density increased, coming closer to the value obtained in the box geometry. 


Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (rnqC) 

Pr,avg (rnqC) 

rmax(mm) 

300 

0 

241.3970 

3.310E-03 

3.847E-06 

5.000949 


1 

239.3932 

1.168E-02 

1.627E-05 

5.009932 


2 

238.6671 

1.774E-02 

1.461E-05 

5.018646 


Table 17. Increased plate radius (20 mm), space charge on. 


Taking it a step further, the plate height was increased to 40 mm and then 60 mm; 
the results are shown in Table 18 and Table 19. As the plate height approaches some 
critical value, the parallel plate geometry comes very close to the box geometry. 
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Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (rneC) 

Pr,avg (rneC) 

rmax(mm) 

300 

0 

259.3666 

3.395E-03 

4.004E-06 

5.000672 


1 

256.5170 

1.195E-02 

1.579E-05 

5.007867 


2 

255.7882 

1.815E-02 

1.331E-05 

5.014946 

Table 18 

Increased plate radius (40 mm), space charge on. 

Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz,avg (rneC) 

Pr,avg (rneC) 

rmax(mm) 

300 

0 

263.0620 

3.421E-03 

3.992E-06 

5.000683 


1 

259.3737 

1.200E-02 

1.514E-05 

5.007839 


2 

258.4446 

1.822E-02 

1.267E-05 

5.015056 


Table 19. Increased plate radius (60 mm), space charge on. 


The results from Tables 14-19 are summarized in Figure 29. 



Plate Radius (mm) 


Figure 29. Variation of current density with plate radius. 


51 















THIS PAGE INTENTIONALLY LEET BLANK 


52 



VII. SIMULATION TIME VARIATIONS 


A. OVERVIEW 

The previous simulations used a simulation time of 0.5 ns. This section 
investigates the impact of simulation time on the current density delivered. 

B. BOX 

If the simulation time is extended, the simple box geometry returns a value 
consistent with the original run time, whether the space charge is on or off. The results 
are summarized in Table 20 (space charge off) and Table 21 (space charge on). 


Time (ns) 

Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 


Time (ns) J 

,ed (kA/m^) 

Screen 

Jdei (kA/m^) 

0.5 

100 

0 

99.83024 


0.5 

300 

0 

299.4907 



1 

99.81605 




1 

299.4482 



2 

99.83352 




2 

299.5006 


100 

0 

99.83331 


1 

300 

0 

299.4999 



1 

99.84319 




1 

299.5296 



2 

99.84277 




2 

299.5283 

2 

100 

0 

99.92712 


2 

300 

0 

299.7814 



1 

99.92187 




1 

299.7656 



2 

99.92384 




2 

299.7715 


Table 20. 

Varying simulation time, box. 

space charge off. 


Time (ns) 

Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 


Time (ns) J 

req (kA/m^) 

Screen 

Jdei(kA/m^) 

0.5 

100 

0 

99.51599 


0.5 

300 

0 

270.2160 



1 

99.31531 




1 

267.0417 



2 

99.28589 




2 

266.0374 

1 

100 

0 

99.70458 


1 

300 

0 

272.0603 



1 

99.60492 




1 

270.7792 



2 

99.59658 




2 

270.9607 

2 

100 

0 

99.86450 


2 

300 

0 

268.5381 



1 

99.80866 




1 

268.4891 



2 

99.80341 




2 

268.7780 


Table 21. Varying simulation time, box, space charge on. 
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C. PARALLEL PLATES 

The parallel plate geometry is affected by the simulation time when the space 
charge is on. As the simulation time increases, the current density decreases. This effect 
is more noticeable when the current density is above the theoretical saturation current 
density. Table 22 shows the results when the space charge is off; Table 23 shows the 
results when the space charge is on. 


Time (ns) 

Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 


Time (ns) 

Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

0.5 

100 

0 

99.83024 


0.5 

300 

0 

299.4907 



1 

99.81605 




1 

299.4482 



2 

99.83352 




2 

299.5006 

1 

100 

0 

99.83331 


1 

300 

0 

299.4999 



1 

99.84319 




1 

299.5296 



2 

99.84277 




2 

299.5283 

2 

100 

0 

99.92712 


2 

300 

0 

299.7814 



1 

99.92187 




1 

299.7656 



2 

99.92384 




2 

299.7715 

Table 22. 

Varying simulation time, parallel plate, space charge off. 


Time (ns) 

Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

0.5 

100 

0 

99.08017 



1 

98.39032 



2 

98.02934 

1 

100 

0 

98.78264 



1 

97.91741 



2 

97.59198 

2 

100 

0 

80.78802 



1 

80.70150 



2 

80.70641 


Time (ns) 

Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

0.5 

300 

0 

187.3599 



1 

184.9992 



2 

185.8662 

1 

300 

0 

136.0386 



1 

136.2075 



2 

136.4122 

2 

300 

0 

85.26074 



1 

85.48097 



2 

85.66611 


Table 23. Varying simulation time, parallel plate, space charge on. 
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VIII. BOX VARIATIONS 


A. OVERVIEW 

At this point, the parallel plate geometry and its variations are set aside. The 
highly variable response of this geometry to increased run times, plate radius, and plate 
thickness preclude setting a reasonable baseline for further study. The following 
simulations examine the impact of varying aspect ratio for the box geometry only. In the 
first set of simulations, the aspect ratio is varied by changing the interelectrode spacing. 
Next, the aspect ratio is held constant while the geometry is scaled up. In all cases, the 
voltage was adjusted to keep the electric field magnitude constant at 1 MV/m. 

B. ASPECT RATIO VARIATIONS 

Changing the aspect ratio by varying the interelectrode spacing also changes the 
theoretical saturation current density. Conversely, changing the aspect ratio by varying 
the cathode radius does not change the theoretical saturation current density. The 
theoretical saturation current changes, but since the current density is determined by 
dividing by the area, beam diameter is removed from the calculation. 

I. AR=200 


The cathode remained at 5 mm while the interelectrode spacing was decreased to 
0.05 mm. 


AR = 


2(^^200 

0.05 mm 


(32) 


The accelerating voltage was adjusted to 50 V to maintain the electric field 
strength at 1 MV/m between the electrodes. The theoretical saturation current density is 
330 kA/rn. 



;r(8.85xlO-'') 12(1.602 x lO'^'^) 


19 N 


(ARf = 


9.109x10 


-31 


(200)^=0.07328 (33) 
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/ = py3/2 ^ (0.07328X50)''' = 25.91 A 


(34) 



25.91 A 
;t( 5 mm)' 



(35) 


Figure 30 shows that the saturation current density delivered by spiffe is 
approximately 380 kA/m'. 



Current Density Requested (kA/m^) 


Figure 30. Box geometry, AR=200. 


2. AR=50 


The cathode remained at 5 mm while the interelectrode spacing was increased to 
0.2 mm. 


0.2 mm 


(36) 


The accelerating voltage was adjusted to 200 V to maintain the electric field 
strength at 1 MV/m between the electrodes. The theoretical saturation current density is 
165 kA/m'. 



;r(8.85xlO-‘') 12(1.602 x lO'*") 


-19 N 


{ARf = 


9.109x10 


-31 


(50)' = 0.004580 


(37) 
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/ = p\/3/2 ^ (0.004580)(200)^'' = 12.95 A 


(38) 


^ 7 12.95 A kA 

J- — - -7 = 165 —7 

A ;t(5 mm) m 


(39) 


Figure 31 shows that the saturation current density delivered by spiffe is 
2 

approximately 180 kA/m . 



Current Density Requested (kA/m^ 


Figure 31. Box geometry, AR=50. 


3. AR=10 


The cathode remained at 5 mm while the interelectrode spacing was increased to 

1 mm. 


AR = 


2(5 mm) ^ 
1 mm 


(40) 


The accelerating voltage was adjusted to 1000 V to maintain the electric field 
strength at 1 MV/m between the electrodes. The theoretical saturation current density is 
73.8 kA/rn. 
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^^0 

9 



;r(8.85xlO-‘^) / 2(1.602 x 10"*^') 
9 i 9.109 xlO''* 


(10)^=0.0001832 
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/ = p\/3/2 ^ (0.0001832X1000)''" = 5.79 A 


(42) 



”^X=73.8 M 

;t(5 mm) m 
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Figure 32 shows that the saturation current density delivered by spiffe is 

2 

approximately 93 kA/m . 
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Figure 32. Box geometry, AR=10. 


4. AR=5 


The cathode remained at 5 mm while the interelectrode spacing was increased to 

2 mm. 


A/? = 


2(5 mm) 
2 mm 


= 5 


(44) 


The accelerating voltage was adjusted to 2000 V to maintain the electric field 
strength at 1 MV/m between the electrodes. The theoretical saturation current density is 
52.2 kA/m". 


^^0 ^(8.85xlO-‘") 2(1.602xl0-‘') 

— (AR) = ---\—fTTP— 
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9 


9.109x10 


4.58x10 


-5 


(45) 
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I = pv3/2 ^ ^ 10-')(2000)^'" = 4.09 A 


(46) 


/ 4.09 A kA 

- =-r = 52.2 — 

A ;t(5 mm) m 


(47) 


Figure 33 shows that the saturation current density delivered by spiffe is 

2 

approximately 67 kA/m . 
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Figure 33. Box geometry, AR=5. 


5. AR=1 


The cathode remained at 5 mm while the interelectrode spacing was increased to 

10 mm. 


AR = 


2(5 mm) 
10 mm 


= l 


(48) 


The accelerating voltage was adjusted to 10,000 V to maintain the electric field 
strength at 1 MV/m between the electrodes. The theoretical saturation current density is 
23.3 kA/ml 



;r(8.85xlO-'') 12(1.602 x 10'*'^) 


-19 > 
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9.109x10 


-31 


(1)' = 1.83x10 
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(49) 
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-^ = 23.3 A 
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Figure 34 shows that the saturation current density delivered by spiffe is 

2 

approximately 40 kA/m . 
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Figure 34. Box geometry, AR=1. 

C. CONSTANT ASPECT RATIO, GEOMETRY SCALED-UP 

In the following simulations, the aspect ratio is maintained at 100 by varying both 
the interelectrode spacing and the cathode radius. The theoretical current density changes 
because the interelectrode spacing has changed. 


1. Scale Up x5 


The cathode radius was increased to 25 mm and the interelectrode spacing was 
increased to 0.5 mm. The grid resolution was adjusted to 0.025 mm. 


AR = 


2(25 mm) 
0.5 mm 


= 100 


(52) 
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The accelerating voltage was adjusted to 500 V to maintain the electric field 
strength at 1 MV/m between the electrodes. The theoretical saturation current density is 
104 kA/m^. 


9 \ m 

^.(8,85 x10-) 2(1.602x10-’) 

9 V 9.109x10'^' 

(53) 


/ = pv^'^ = (1.83 X 10‘")(500)^'' = 205 A 

(54) 


/ 205 A M 

/ = —=-7 = 104 — 

A .;t( 25 mm) m 

(55) 


Figure 35 shows that the saturation current density delivered by spiffe is 

2 

approximately 118 kA/m . 
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Figure 35. Box geometry, magnification x5 

2. Scale Up xlO 

The cathode radius was increased to 50 mm and the interelectrode spacing was 
increased to 1 mm. The grid resolution was adjusted to 0.05 mm. 
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^ „ 2(50 mm) 

AR = — -^ 

1 mm 


= 100 


The accelerating voltage was adjusted to 1000 V to maintain the electric field 
strength at 1 MV/m between the electrodes. The theoretical saturation current density is 
73.7 kA/ml 


9 \ m 

^.(8.85x10-) 2(1.602x10-^ 

9 V 9.109x10'^' 

(57) 


/ = pv 3/2 ^ g3 ^ io-')(1000)''' = 579 A 

(58) 


=73.7*1 

A ;t(50 mm) m 

(59) 


Figure 36 shows that the saturation current density delivered by spiffe is 

2 

approximately 78 kA/m . 
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Figure 36. Box geometry, magnification xlO 

3. Scale Up x50 

The cathode radius was increased to 250 mm and the interelectrode spacing was 
increased to 5 mm. The grid resolution was adjusted to 0.25 mm. 
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The accelerating voltage was adjusted to 5000 V to maintain the electric field 
strength at 1 MV/m between the electrodes. The theoretical saturation current density is 
33.0kA/ml 


9 \ m 

^.(8.85x10-) 2(1.602x10-^ 

9 V 9.109x10'^' 

(61) 


/ = pv 3/2 ^ g3 ^ io-")(5000)''' = 6477 A 
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Figure 37 shows that the saturation current density delivered by spiffe is 

2 

approximately 25 kA/m . 



I 

e 

9 

o 

9 

b 

a 



Current DensItY Requested (kA/m^) 


Figure 37. Box geometry, magnification x50 

4. Scale Up xlOO 


The cathode radius was increased to 500 mm and the interelectrode spacing was 
increased to 10 mm. The grid resolution was adjusted to 0.25 mm. 
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AR = 


2(500 mm) 
10 mm 


The accelerating voltage was adjusted to 10,000 V to maintain the electric field 
strength at 1 MV/m between the electrodes. The theoretical saturation current density is 
23.3 kA/ml 


p _ ^^0 


—(ARf 
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9.109x10 
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I = pv 3/2 ^ ^ io-^)(iooOO)"'" = 18320 A 
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(65) 

( 66 ) 


/_ 18320A _333M 
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Figure 38 shows that the saturation current density delivered by spiffe is 

2 

approximately 14.5 kA/m . 
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Figure 38. Box geometry, magnification xlOO. 
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IX. ION SIMULATION 


A. PROGRAM SET-UP 


Spiffe provides one option for altering the program to simulate ions. This is 
accomplished through use of the “stiffness” command. The stiffness is particle mass, 
given in terms of electron mass. For Argon, the stiffness is 72819, meaning that Argon is 
72819 times as massive as an electron. 


39.948 amu 


^1.6605x10 


amu 


J K9.1094 X IQ-"'kgJ 


= 72819e" 


( 68 ) 


There is no provision for changing the charge of the particle, so the problem must 
be set up in reverse. It is also helpful to increase the accelerating voltage. Since the 
particle is now significantly more massive, the acceleration will decrease for the same 
accelerating voltage (force) applied. This would require a longer simulation time for the 
particles to traverse the interelectrode spacing. In order to keep the simulation time at the 
same value used in the baseline simulations (0.5 ns), the potential was increased by the 
same order of magnitude as the stiffness. The potential of the original box geometry was 
increased from 100 V to 7 MV. If the potential had not been increased, the ion 
simulation would have to run 10^ times longer to release the same number of 
macroparticles as the electron simulation. This begins to approach the maximum 
simulation time, limited by 32-bit signed integer memory. 


B. RESULTS 


With the accelerating voltage at 7 MV and AR=100, the theoretical saturation 

2 

current density is 16 GA/m . 


P = = ^(8-85x10-) 2(1.602 xlO-J) ^ ^ ^ 

9 \ m 9 V 6.633x10-'' 


(69) 


/ = pv^'^ = (6.79 X10"^ )(7 X 10'’ = 1.26 MA 


(70) 
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Figure 39 shows that the saturation current density delivered by spiffe is 
approximately 18.5 GA/m^. 
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Figure 39. Ion simulation, box geometry, AR=100. 


















X. ANALYSIS 


A. ONE-DIMENSIONAL APPROXIMATION 

Spiffe accurately models the theoretical saturation current density in a one¬ 
dimensional approximation. The current density delivered by the simulation is 
approximately 15% higher than the theoretical value. 

The slightly higher value of current density is expected due to beam spreading. 
The Child-Langmuir law applies to infinite parallel plates. In a beam of finite width, the 
space charge interaction between particles will cause beam spreading. The high aspect 
ratio of the box approximates a one-dimensional situation, but because the beam is finite, 
it will still experience spreading. 

The electric field from a point charge is a function of 1/r . In a two-dimensional 
scenario with a beam of finite width, interactions between the ions cause the beam to 
spread radially. As the beam spreads (i.e., r increases), the effect of the ions’ electric 
field at the source is diminished. Since the electric field from the ions counteracts the 
external electric field, a decrease in the ion electric field allows an increase in the external 
electric field. The increased external electric field allows the source to emit a current 
density that is higher than the theoretical limit given by the Child-Langmuir law. 

Beam spreading in the parallel plates is greater, which implies a higher saturation 
current density, but this was not the result. The average longitudinal momentum in the 
parallel plate geometry is lower than that seen in the box geometry. This indicates that 
although the interelectrode spacing and applied potentials are the same, the parallel plate 
geometry is not accelerating particles as much as the box geometry. 

In both the box and parallel plate geometries, with the space charge off, spiff e 
returns the expected current density and electric field distribution. With the space charge 
on, the geometries vary significantly. The box geometry returns a higher saturation 
current density than expected, which is explained by beam spreading. In contrast, the 
parallel plate geometry returns a much lower current density than expected. The electric 
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field plots provide some insight into this difference. Before particles are released, the 
equipotential lines of the box geometry are uniformly spaced. The equipotential lines of 
the parallel plate geometry are concentrated at the cathode and anode. The magnitude of 
the electric field at the centerline (r=0) of the box ranges from 4x10^ V/m to 1.4x10^ 
V/m. The range in the parallel plate geometry begins at 4x10^ V/m but only increases to 
1x10^ V/m. This is perhaps why the particles are not being accelerated as in the box 
geometry. The cause for these differences is likely due to the difference in boundary 
conditions for the cavity border and the physical surface inside the cavity. Both 
boundary conditions require E^=0, but the program interprets them differently when the 
space charge calculation is involved. 

Removal of one plate from the parallel plate geometry confirms that it is the 
boundary conditions causing the differences between the box and parallel plate 
geometries. The variability of current with plate thickness and radius is another indicator 
that the boundary conditions are the cause for variations between the geometries. 

As the plate radius increases, it approaches the values provided by the box 
geometry. The excessive radius needed to reach this situation is unnecessary when the 
simpler box geometry can accomplish the same result. Variations between the box and 
parallel plate geometries illustrate that the cathode should be coincide with a cavity 
border. 

B. BOX VARIATIONS 

In the one-dimensional situation, the current density was approximately 15% 
higher than the theoretical limit. As the aspect ratio is decreased, the current density 
deviates further from the theoretical limit (Table 24, Figure 40). This is not necessarily 
an error in the modeling, but a result of the fact that beam spreading becomes significant 
at lower aspect ratios. The practical aspect ratio limit for a beam is on the order of unity, 
as noted in [2]. 
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AR 


200 

100 

50 

10 

5 

1 

Theoretical (kA/m ) 

329.89 

233.26 

164.94 

73.76 

52.16 

23.33 

■J delivered (kA/m ) 

381.25 

268.64 

181.54 

93.54 

66.19 

39.50 

■Jdel/Th 

1.16 

1.15 

1.10 

1.27 

1.27 

1.69 


Table 24. Theoretical vs. delivered current density, varying AR. 
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Figure 40. Variation in current density delivered/theoretical for varying aspect ratios. 


Spiffe is best used for small-scale geometries. Table 25 compares the current 
density for the magnified geometries. With the aspect ratio at 100 in all cases, spiffe 
should have returned a current density 15% greater than the theoretical limit. As the 
scaling of the geometry increased, the current density decreased (Table 25, Figure 41). 


AR=100 

Magnification 

1 

5 

10 

50 

100 

Theoretical (kA/m ) 

233.3 

104.3 

73.76 

32.99 

23.33 

^delivered (kA/m ) 

268.6 

117.1 

77.65 

24.37 

14.37 

Tel/Th 

1.15 

1.12 

1.05 

0.74 

0.62 


Table 25. Theoretical vs. delivered current density, AR=100, geometry magnified. 
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Figure 41. Variation in current density delivered/theoretical for varying 

magnific ations, AR= 100. 


C. ION SIMULATION 

With AR=100 and the simulation adjusted to simulate ions, the saturation current 
density is approximately 15% higher than the theoretical value. This is the same ratio 
that was seen for the electron simulation. The voltage was increased significantly to 
reduce simulation time. Voltage breakdown is a concern, but not at this stage. 
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XI. CONCLUSION AND RECOMMENDATIONS 


A. CONCLUSION 

Spif f e differs from current steady-state ion numerical modeling schemes in that 
it solves Maxwell's equations and performs a Poisson correction. As the program was 
intended to simulate electrons, many quantities are normalized to the mass and charge of 
an electron. 

The code is able to simulate ions by adjusting the mass of the particle, but it 
provides no provision for changing the charge. This limits the simulation particles to 
singly-charged ions. This is satisfactory for field-ionization, which produces only singly- 
charged ions. The sign of the charge remains negative, so the user must be prepared to 
set-up the problem in reverse. 

Spif fe accurately models a theoretical one-dimensional situation. Finite beam 
spreading causes the program to output a current density approximately 15% higher than 
the theoretical limit. As the geometry deviates from one-dimensional, the error increases 
due to radial beam spreading. This is not necessarily an error in the modeling, but a 
result of the fact that beam spreading becomes significant at lower aspect ratios. 

A CNT field-ionized ion thruster model should be capable of modeling the 
electrode spacing and applied potentials to both interact with the CNTs to create ions and 
also to accelerate those ions. Spiffe is not capable of doing this for two reasons. The first 
is that spiffe cannot solve the Fowler-Nordheim equation, making it unable to model 
ionization of the propellant gas. Second, the user would be required to program the 
geometry of every nanotube, and the grid would have to be small enough to resolve those 
nanotubes and fill the entire accelerating region. The time step required for stability 
would be small and the simulation time would be long, both contributing to impractical 
run times and CPU memory limitations. 

Spiffe is a useful tool for the early modeling stages of a CNT field-ionized ion 
thruster. It can be used to predict the saturation current density for the setup described in 
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[8]. Additionally, it can be used to model the ideal eleetrode spaeing for a single beamlet 
with a simplified planar ion souree. Sinee it is strietly an axi-symmetrie eode, spiffe 
is not able to model a multi-aperture grid. As the CNT field-ionized ion thruster is 
further developed, a more robust model should be pursued. 

B. FUTURE WORK 

Future work ineludes using spif f e to prediet the current producing capability of 
the setup deseribed in [8]. 

Spif fe may also be used for determining the thrust and perveanee for a single 
beamlet in order to eonduct feasibility studies for CNT field-ionized thrusters and micro- 
and nano-seale satellites. 

Finally, future work must investigate alternate modeling options for full-seale 
numerieal modeling of CNT field-ionized thrusters for miero- and nano-seale satellites. 
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APPENDIX A. SUPPLEMENTAL USER’S GUIDE 


A. USER’S GUIDE 

Reference [14] is the user’s guide published by Argonne National Laboratory 
(ANL)/Advanced Photon Source (APS). It is available for download from the ANL 
website. This Appendix is intended to supplement the official user’s guide to provide 
additional guidance for first-time users. 

B. SDDS TOOLKIT 

SDDS Toolkit is a set of commands that allows plotting and processing of the 
data obtained from spiffe runs. The tool ki t and its user’s guide are available for 
download from the ANL/APS website. 

C. DOWNLOAD AND INSTALL 

1. Download Address. 

The ANL/APS software download page is located at: 

http://www.aps.anl.gov/Accelerator_Systems_Division/Accelerator_Operations_P 

hysics/software.shtml 

2. Spiffe 

Scroll down to the spiffe section and download the appropriate file for your 
operating system. The simplest method is to use the self-installing file (Spiffe.exe). 

3. SDDS Toolkit 

Scroll down to SDDS Binaries and download the appropriate file for your 

operating system. The simplest method is to use the self-installing file (SDDS 
ToolKit.exe). 
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4. 


Verify Installation 


Once the spiff e installation is complete, open a command prompt window and 
type “spiffe” at the prompt. If installed correctly, the prompt will return the following 
lines: 

Link date: Sep 6 2011 08:52:07 (or date of latest revision) 

error: no input file listed 

usage: spiffe inputfile 

Program by M. Borland (This is version 4.1, August 2011.) 

D. PROGRAM STRUCTURE 

Spiffe files may be written in any text editing program. A spiffe run 
consists of at least two files - the geometry file and the main run file. The geometry file 
is an input to the run file. The run file may also output data to several files, which are 
specified by the user. The run file must have a “.spiffe” extension. The geometry file 
and all output files specified by the user may be given any extension. Common practice 
is to use something that indicates the type of data contained in the file. For example, 
“.geo” for the geometry file and “.snap” for the beam snapshot data. The recommended 
sequence of commands for the run file is given in [14]. 

E. GEOMETRY FILE 

I. Overview 

The geometry file can be given any name and any extension. The geometry file 
defines the geometry using a series of points. Each point may be assigned a potential. If 
a potential is not specified, the default is 0 V. Any geometry that can be defined with line 
segments and arcs can be defined with the geometry file. The most important thing to 
remember when defining a geometry file is that spiffe is an axisymmetric code. The 
geometry file uses x-y coordinates to define each point. When translated by the run file, 
v-coordinates correspond to the z-axis and y-coordinates correspond to radius. 

The geometry file may be written using any unit of length, but meters is preferred 
for simplicity. If another unit is used, the spiffe run file provides a command line 
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(lr_factor) to convert units to meters. 

An optional command that is not included in [14] is ramp_potential. This 
command allows you to ramp the assigned potential along the surface to prevent 
discontinuities. In the simple geometries used in this report, ramp_potential made 
little difference in the results. It may become important when defining boundaries inside 
the cavity. 

2. Example - Box Geometry 

The structure of the geometry file is described in Section 4.2 of [14]. This 
example illustrates how to create the box geometry used in this report. 

First, open a new text editor file and name it box.geo. Type the following code 
into the text file: 

&po x=0.0, y=0, potential=0, &end 

&po x=0.0, y=0.01, potential=0, ramp_potential=l, &end 

&po x=0.0001, y=0.01, potential=100, ramp_potential=l, &end 

&po x=0.0001, y=0.0, potential=100, &end 

&po x=0.0, y=0.0, potential=0, &end 

Each line represents a new point. Note that &po is an acceptable abbreviation 

for &point. Five points are defined because it is points that are defined, not line 

segments. If the first point was not restated, the fourth line segment would not be drawn. 

Ramp_potential is not used at y=0 since this is the centerline (r=0) of the 

axisymmetric beam. 

F. RUN FILE 

1. Overview 

The run file calls on the geometry file, draws the geometry, and defines all aspects 
of the simulation (e.g., cathode radius, run time, integration step size, grid size, etc.) 

The first step in creating a spiffe run file is to establish the geometry for the 
simulation. This is done with the geometry file and the def ine_geometry command, 
described in section 4.1 of [14] and discussed below. 
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2 . 


Define the geometry 


a. Cavity 

The simulation takes place in a cavity. In the case of the box geometry, 
the simulation cavity is the box. In the case of the parallel plate geometry, the cavity is 
the region surrounding the plates. The simulation cavity is defined by specifying the 
minimum and maximum longitudinal coordinates (zmin, zmax) and the maximum 
radial coordinate (rmax). The minimum radial coordinate is always r=0 since spiffe 
is an axisymmetric code. 

b. Grid 

The grid is defined by defining the cavity (zmin, zmax, rmax) and 
specifying the number of grid lines (nz, nr). Grid lines are equally spaced between the 
minimum and maximum values; the shape is a simple quadrilateral. Radial and 
longitudinal grid resolution should be kept within a factor of two for best results. 

A knowledge of the grid layout is essential for drawing geometries and 
assigning the proper number of grid lines. It is also helpful for interpretation of electric 
field plots and understanding how the space charge is calculated. 

Field quantities are calculated using four interleaving grids. The grid and 
position-index correspondences are described in [12]. Figure 42 illustrates the 
relationship between the four interleaving grids. Black lines indicate the grid that is 
specified by the user with nz and nr, which is also the grid for calculating p and O. The 
longitudinal components of the electric field (Ez) and the current density (Jz) correspond 
to the red grid (horizontal lines overlap black). The radial components of the electric 
field (Er) and current density (Jr) correspond to the blue grid (vertical lines overlap 
black). Einally, the grid for the magnetic field (B,|,) overlaps the red and blue lines. 
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Figure 42. Interleaving grid. 


In the box geometry, the only dimensions of eoneern are the length (0.1 
mm) and radius (10 mm) of the simulation region. The longitudinal and radial grid 
spacing are defined in [12] as: 


Az 


z — z ■ 

^max ^mm 


Ar = 


r — r ■ 

max _ mm 

n,-l 


For a grid resolution of Az = Ar = 0.05 mm. 


0.1 mm 
0.05 mm 


+ 1 = 3 


nr - 


10 mm 
0.05 mm 


+ 1 = 201 


(72) 

(73) 


(74) 

(75) 
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1 ^ 

With this grid size, spiffe requires a time step less than 1.179327x10' sec for stable 
integration. 


Table 26 compares spiffe output for various grid sizes and current 
density. The grid size has minimal impact when the requested current density is less than 
the saturation current density. The saturation current for the box geometry is 
approximately 233 kA/m . With the largest grid resolution (fewest number of grid lines), 
the simulation returns a value well beyond the space charge limit. As the resolution 
improves (grid lines increase), spiffe provides a better approximation of the space 
charge limit. 


Resolution 

50 pm 

25 pm 

10 pm 

5 pm 

z grid lines (nz) 

3 


5 


11 


21 

r grid lines (nr) 

201 

401 

1001 

2001 

Maximum time step 

1.18E-13 

5.90E-14 

2.36E-14 

1.18E-14 

Time step used 

lE-13 

5E-14 

2E-14 

lE-14 

Number of steps 

5000 

10000 

25000 

50000 

Jreq=100 kA/m^ 

-Idel ^max 

-Idel 

tmax 

-Idel 

tmax 

-Idel 

tmax 

Screen 0 

99.855 4.999998 

99.288 

4.999986 

99.674 

5.000044 

99.517 

5.000102 

Screen 1 

99.908 5.000202 

99.086 

5.000581 

99.508 

5.001042 

99.312 

5.001586 

Screen 2 

99.950 5.000606 

99.080 

5.001537 

99.429 

5.002419 

99.283 

5.003382 

Jreq=300 kA/m^ 

-Idel ^max 

-Idel 

tmax 

-Idel 

tmax 

-Idel 

tmax 

Screen 0 

299.04 5.000016 

296.32 

5.000058 

291.53 

5.000261 

270.14 

5.000620 

Screen 1 

298.42 5.000761 

293.55 

5.002499 

288.34 

5.004768 

266.93 

5.007441 

Screen 2 

298.27 5.002157 

293.47 

5.005997 

287.49 

5.010420 

265.92 

5.014438 

Jreq=500 kA/m^ 

-Idel ^max 

-Idel 

tmax 

-Idel 

tmax 

-Idel 

tmax 

Screen 0 

497.64 5.000045 

405.78 

5.000223 

307.71 

5.000858 

271.10 

5.001814 

Screen 1 

494.12 5.001426 

402.27 

5.007025 

300.65 

5.012518 

268.08 

5.013434 

Screen 2 

490.98 5.004336 

398.27 

5.015544 

299.61 

5.022965 

266.73 

5.023683 


Table 26. Impact of grid size on simulation accuracy. 


There is a tradeoff between simulation time and grid resolution. As the 
grid lines increase, the time step required for stability decreases, which increases the run 
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time of the simulation. In the simplest case (3x21), the simulation took less than 10 
seconds. With the grid size used throughout this report, simulations averaged 20 to 40 
minutes. 


In the parallel plate geometry, the smallest dimension is the plate thickness 
at 0.01 mm (10 pm). The grid must be able to distinguish plate boundaries, so it must 
have a minimum resolution of 10 pm. The grid must also provide enough resolution to 
accurately model the particle and field interactions between the plates. The gap between 
the plates is 0.1 mm; a resolution of 10 pm would place ten grid lines between the plates. 

The cavity surrounding the plates is 0.14 mm wide and 10 mm tall. Using 
a resolution of 10 pm (needed to distinguish plate thickness) and Equations (72) and (73), 
the number of longitudinal and radial grid lines is calculated by: 


0.14 mm 

nz^ - 

0.01 mm 


+ 1 = 15 


(76) 


nr - 


10 mm 
0.01 mm 


+ 1 = 1001 


(77) 


When using this grid resolution, spiffe reports that the time step must 
be less than 2.540089x10'^^ for stable integration. 


c. Geometry Output Files 

Spiffe writes two output files regarding the geometry. The first file is 
the geometry as it is specified by the user (boundary_output). The second file is the 
geometry as it is interpreted by spiffe (discrete_boundary_output). These 
files are helpful for verifying that the geometry was programmed or interpreted correctly. 

It is essential to choose grid lines that coincide with boundaries. This can 
be verified by comparing the ideal boundary (user-defined) to the discrete boundary 
(spiff e-interpreted). For example, in the parallel plate model, the grid lines match up 
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with the plate boundaries (Figure 43). If the grid lines were chosen poorly, they would 
not line up with the plate boundaries and spiffe would not “see” the plate correctly. 
This is illustrated in Figure 44. 

- 1 - 1 - 1 - 

0,01 0 . 

0.008 _ 

^ 0.006 _ 

^ 0,004 . 

0,002 . 

0.000 . 

0 5.0x1 0-5 1 .0x1 0-4 

z (m) 

Figure 43. Geometry interpretation with correct number of grid lines. 

0.01 0 

0.008 
^ 0.006 
0.004 

0.002 

0.000 


Figure 44. Geometry interpretation with incorrect number of grid lines. 

SDDS Toolkit is needed to produce the plots in Figure 43 and Figure 44. 
If the ideal boundary file is named duo.bnd and the discrete boundary file is named 
duo.dbnd, the following command will produce the comparison plot. 



0 5.0x1 0-5 1 .0x1 0-4 

z (m) 
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sddsplot -col=z,r -graph=line,type=2 -split=page duo.bnd - 
col=z,r -graph=line,type=3 -split=page duo.dbnd 

Further details of the SDDS plot tools can be found in [15]. 

d. Boundary Conditions 

Boundary conditions for the cavity are specified in the 
define_geometry section of the code. Upper, lower, left, and right boundary 
conditions apply to the boundaries of the cavity. 

Boundary conditions for surfaces inside the cavity (e.g., plates) are 
assumed to be metal surfaces. In this case, £^=0 for vertical surface and Er=0 for 
horizontal surfaces. In other words, there are no electric field lines normal to the metal 
surfaces. 


e. Radial and Longitudinal Interpolation 

Radial_interpolation and longitudinal_interpolation 

should be set to “1” to indicate that field components should be interpolated when fields 
are applied to particles. This is necessary to prevent an abrupt change as particles move 
from one grid square to the next [14]. 

/. Example Code 

The code used for the box geometry is included here. The code is written 
by simply opening a new text file and copying the text shown. The unit of measurement 
for everything specified in a spif fe run file is meters. The command zr_f actor is a 
scaling factor that is used to convert the geometry file (box.geo) into meters. Since 
box.geo is already written in meters, the command zr_factor can be set to “1” or 
simply omitted from the code. Any items not included in the code are set to their default 
value. Additional code options and their default values are listed in Section 4.1 of [14]. 


&define_geometry 
nz=21, 
zmin=0, 
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zmax=0.0001, 

nr=2000, 

rmax=0.01, 

boundary="box.geo", 

radial_interpolation=l, 

longitudinal_interpolation=l, 

boundary_output="%s.bnd", 

discrete_boundary_output="%s.dbnd", 

lower="Dirichlet", 

upper="Neumann", 

right="Neumann", 

left="Neumann", 

Send 


If the name of the spiffe run file is box.spijfe, then the notation %s will 
send the ideal boundary points to box. bnd and the discrete boundary points to box. dbnd. 
The parallel plate geometry was created with the following code. 


! Define the cavity 


&po 

x=- 

-0.00002, 

II 

o 

o 

potential=0. 

Send 

&po 

x=- 

-0.00002, 

y=0.011 

, potential=0. 

Send 

&po 

x= 

0.00012, 

y=0.011 

, potential=0. 

Send 

&po 

x= 

0.00012, 

II 

O 

o 

potential=0. 

ramp_potential=l, 

&po 

x= 

0.00011, 

II 

o 

o 

potential=100. 

ramp_potential=l, 

&po 

x= 

0.00010, 

II 

o 

o 

potential=100. Send 

&po 

x= 

o 

o 

II 

o 

o 

potential=0. 

Send 

&po 

x=- 

-0.00002, 

II 

o 

o 

potential=0. 

Send 


Send 

Send 


IDefine the left plate 


&po 

nt=3, x=0.0. 

II 

O 

o 

potential=0. 

Send 

&po 

X 

II 

o 

o 

t—1 

o 

o 

II 

>1 

potential=0. 

Send 

&po 

x=-0.00001, 

\—1 
o 

o 

II 

>1 

potential=0. 

Send 

&po 

x=-0.00001, 

II 

o 

o 

potential=0. 

Send 

&po 

X 

II 

o 

o 

II 

o 

o 

potential=0. 

Send 

IDefine the right plate 



&po 

nt=3, x=0.0001. 

II 

o 

o 

potential=100. 

Send 

&po 

x=0.0001, 

t—1 

o 

o 

II 

>1 

potential=100. 

Send 

&po 

x=0.00011, 

t—1 

o 

o 

II 

>1 

potential=100. 

Send 

&po 

x=0.00011, 

o 

o 

II 

potential=100. 

Send 

&po 

x=0.0001, 

II 

o 

o 

potential=100. 

Send 


3. Define the Cathode 

This section of the code is where the user specifies the location of the electron 
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source, the current density requested, and the properties of the particles when they 
initially appear in the simulation (e.g., initial momentum). 

a. Position 

The physical cathode is defined in the geometry file, but there is no 
requirement to place the particles on the physical surface. Particles may be loaded from 
any location that can be described by radius and longitudinal position. Placing the 
cathode slightly off the physical surface theoretically ensures that no particles are 
artificially lost to the boundary. In reality, the cathode position (z_posit ion=0 versus 
z_posit ion=le-8) had minimal impact on the current delivered. As the size of the 
cathode increases, this may become more significant. 

b. Current Density 

The user can request any current density. If the value requested is greater 
than the saturation current, spiffe will not exceed the saturation current (within 
numerical error and considering finite radius beam spreading). Particles are generated 
randomly from the location specified. 

c. Macroparticles 


Spiffe is a particle-in-cell (PIC) code. PIC codes are characterized by 
the use of macroparticles, which represent several physical particles. The user may 
specify either electrons_per_macroparticle or number_per_step. 
Electrons_per_macroparticle tells spiffe how many electrons to assign to 
each macroparticle. For example, if the user specifies 500 electrons per macroparticle 
with a time step of 1x10'^"^ sec, a current density of 5000 A/m^, and a cathode radius of 5 
mm, spiffe will calculate the appropriate number of macroparticles (mp) to release per 
time step. 


50004 - 


HI) C 

A-sIC 1.602xl0-'‘'C 


IxlO" 


= 49 ^ 
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If instead the user speeifies the number of maeroparticles to release per 
time step, spiffe will ealeulate the number of eleetrons that are represented by eaeh 
maeropartiele. 


5000 A - e- 

A-sic 1.602x10'^C step 8 mp mp 


(79) 


The ehoiee is entirely up to the user. Fewer eleetrons per maeropartiele 
will improve aeeuraey slightly but drastieally inerease eomputation time and memory. 

d. Start/Stop Time 

These times apply only to the release of particles. The simulation time is 
specified in a different section and does not need to match cathode start and stop times. 
Generally, it is recommended to start the cathode shortly after the simulation starts. This 
allows the initial external electric field to be established before releasing particles into the 
cavity. 


e. Initial Momentum 

The user selects the initial momentum, both longitudinal (init_pz) and 
radial (init_pr), to give the particles once they are randomly generated. This option 
can have a significant impact on the results and must be chosen carefully to properly 
represent the physical scenario. Initial momentum was set to zero for all simulations in 
this report. 

Momentum in spiff e is dimensionless, normalized to electron mass and 
the speed of light (ruec). If the initial momentum is specified as 1 xlO'^, this is equivalent 
to: 

1X 10^'m^c = (1X 10'')(9.109 x 10'‘/tg)(^2.998 — j = 2.731 x 10'®kg-m/s (80) 
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/. Stiffness 


Spiff e allows the user to alter the mass of the particle by adjusting 
beam stiffness. Stiffness refers to the particle mass in terms of electron mass. When 
simulating ions, the stiffness is increased to the appropriate value, but the charge will still 
be negative. The problem is simply set-up in reverse. 

g. Current Density Profile 

Spiffe provides an option to load a current density profile, allowing the 
user to gradually load a current density (e.g., a half-gaussian or linear ramp-up to the 
desired current density). This may be helpful for suppressing numerical oscillations 
while spiff e is converging to the solution. Use of this function requires caution when 
calculating the current density delivered. The screens that collect beam data must start 
after the ramp-up is complete. If the screens collect beam data during the ramp-up, then 
the calculated current density will appear low. 

Minimal work was done with this function because the screens failed to 
start when specified. The current density was ramped linearly from t=0 to t=lxlO'^^ sec 
and then held constant from 1x10'^^ sec to 5x10'^^ sec. The screens were coded to start 
collecting data at 2x10'^^ sec. The calculated current density from the screen was 
approximately 233 kA/m , which happens to be the theoretical saturation current density. 
Closer examination revealed that the screen data included the linear ramp-up, so values 
from J=0 to J=300 kA/m were also included in the calculation. When the raw data was 
extracted from the simulation, the first data point in the first screen corresponded to 1x1 O' 
sec, when it should have been on or near 2x10'^^ sec. After eliminating the undesirable 
points {t < 2x10' sec), the calculated current density was approximately 270 kA/m . 

h. Example Code 

The code used for both the box and parallel plate geometries is shown 
here. Current density was varied manually for separate trials, although this can be done 
automatically with additional programming tools. Number_per_step and 
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start_time were constant for all trials. Time is given in seconds; positions are given 
in meters. This section of code should be added below the def ine_geometry section 
of the text file. 

&define_cathode 

z_position=0, 

outer_radius=5e-3, 

current_density=300e3, 

start_time=le-l1, 

stop_time=2e-9, 

number_per_step=8, 

initial_pz=0. 

Send 

4. Poisson Correction 

The field propagation routine uses two of Maxwell’s equations: Ampere’s law 
and Faraday’s law. It does not use the divergence equation. Accumulated errors in the 
field propagation routine lead to non-zero divergence. Spiffe corrects for this by 
solving the Poisson equation at user-defined intervals during the simulation. Varying the 
options in this section of the code had minimal impact on the simulation results. 

a. Correction Interval 

The settings used in the Poisson correction section of the code have 
minimal impact on simulation time and current density output. The correction interval 
used throughout this report was 32. If the correction interval is increased to 320, the 
simulation finishes about two minutes faster. The current density delivered is similar - 
270.7 kA/m^ versus 271.3 kA/m^ for 32 and 320, respectively. 

b. Accuracy 

Poisson correction accuracy is the fractional accuracy of the Poisson 
solutions. Its impact on the simulation is minimal. With the fractional accuracy set to 
1 xlO'"^, spiffe returned a current density of 270.743 kA/m^. Increasing the fractional 
accuracy by two orders of magnitude caused a negligible change in the output, 270.740 
kA/ml 
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c. Iterations 

The default number of iterations is 1000. If spif f e does not converge to 
a solution within the defined number of iterations, the program will provide a warning. 

Warning: maximum number of iterations exceeded in 

solve_poisson_cyl—residual was 3.952478e+001. 

The simulations used throughout this report left the maximum number of 
iterations at the default value. The impact of this option on the program results is 
minimal. For comparison, if the number of iterations was limited to 100, spiffe 

o 

provided approximately 80 warnings, each with a residual on the order of 10 . With 300 

2 2 

kA/m requested, spiffe delivered 270.809 kA/m . With the default number of 

2 

iterations, spiffe delivered 270.216 kA/m . 

d. Example Code 

The code below was used for all simulations in this report, regardless of 
geometry, run time, and current requested. 

&poisson_correction 

step_interval=32, 
accuracy=le-4, 
verbosity=0, 

error_charge_threshold=le-15. 

Send 

5. Screens 

The define_screen command is described in Section 4.13 of [14]. This 
option is necessary to determine current density output. With the simple geometries in 
this report, none of the screens observed backward particles. In more complicated 
geometries, it may be necessary to record both forward and backward particles. Screens 
may be placed individually using separate commands, or one command may be used to 
place several screens at regular intervals. Position and start time are in meters and 
seconds, respectively. The data is saved to the file that is specified by the user. In the 
example below, the file will be saved to box.scO. Each def ine_screen command 
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should be given a different extension to avoid overwriting data (e.g., “.scO”, “.scl”, 
“.sc2”). 

&define_screen 

filename="%s.scO", 
z_position=le-5, 
direction=forward. 

Send 


The screen output file (e.g., box.scO) contains information on every macroparticle 
that passes through the screen. 

• Radius (r) in meters; the radius of the particle when it crosses the screen. 

• Momentum (pz, pr, pphi) in terms of mgC, the momentum of the particle 
when it crosses the screen. 

• Time (t) in seconds; time that the particle crosses the screen. 

• Charge (q) in Coulombs; will be a multiple of electron charge, 
corresponding to the number of electrons per macroparticle. 

• Initial radius (rO); radius that the particle was initially generated. 

• Initial time (tO); time that the particle was initially generated. 

• particlelD; each macroparticle is assigned an ID number when it is 
generated. 

• Minimum and maximum positions (rmin, rmax, zmin, zmax); minimum 
and maximum position of the particle from initial generation until it 
crosses the screen. 

6. Snapshots 

Beam snapshots are described in Section 4.15 of [14]. Beam snaphots are 
necessary for viewing particle trajectories. They are also helpful for troubleshooting 
unexpected results or problematic code. 

In the following code, the beam data is saved to box.snap. The output file may be 
given any extension.. 

&define_snapshots 

filename="box.snap" 
time_interval=5e-l1, 
start_time=0. 

Send 
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The snapshot output file contains data on every macroparticle for each snapshot in 

time: 

• Position (z, r) in meters. 

• Momentum (pz, pr, pphi) in terms of rUeC. 

• Charge (q) in Coulombs; will be a multiple of electron charge, 
corresponding to the number of electrons per macroparticle. 

• Initial radial position (rO); the radius at which the particle was initially 
generated (random). 

• Initial axial position (zO); the longitude at which the particle was initially 
generated (should correspond to z_position from 
def ine_cathode). 

• Initial time (tO); time that the particle was initially generated. 

• particlelD; each macroparticle is assigned an ID number when it is 
generated. 

7. Field Output Options 

a. Field Output. 

Spiff e provides two methods of examining the electric field. The first, 
def ine_f ield_output, maps the electric and magnetic field of the entire simulation 
area. This is described in Section 4.16 of [14]. The output file may be given any name 
and any extension. Individual data will rarely be pulled from this output file; field plots 
are more useful. Field plots are obtained using the sddscontour command. 

b. Field Sampling 

The second method of examining the field is through 
de f i ne_f i e 1 d_s amp ling. This command provides several options for examining a 
field component at a given position. The output file may be given any name and 
extension desired. Standard practice was to use an extension that described the 
component being sampled and the direction along which it was sampled. For example. 
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when sampling along z at the centerline, the extension was “.Ez_z”. The option used 
throughout this report was to sample along the centerline. This was accomplished 
with the following code: 


&define_field_sampling 
filename="%s.Ez_z" 
component=Ez, 
min_coord=0, 

max_coord=0.0001 (0.000105) 
position=0, 
direction="z", 
time_interval=5e-ll, 
start_time=0, 

&end 


When plotted, the E^ data points will correspond to the E^ grid rather than the 
user-defined grid. If the grid lines (defined by nz) are at z=0 and z=5 mm, then the first 
E^ grid point will lie at z=2.5 mm. In the case of the parallel plate geometry, 
max_coord was set to 0.105 mm to capture the final E^ grid point. 

The field sampling command also provides an option for determining the average 
value of the field along a line. This option was not used for this report. This function did 
not work when requesting a time average along r. 

8. Integration 

a. Time Step 

The simulation step size (dt_integration) is determined by the 
number of grid lines. If the step size is too large for the grid, spiffe will provide a 
warning. 

A small time step will preclude longer simulation times due to memory 

g 1 

limitations. The total number of steps is limited to 2 steps (32-bit signed integer). If 
the time step is IxlO'^"^ seconds, then the simulation is limited to approximately 21.5 ps. 
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231 steps X- ^—- = 2.15 X 10 ^s 

step 


( 81 ) 


Smaller step sizes will increase the run time. More gridlines require a 
smaller step size. There is a tradeoff between accuracy, memory, and time. 

b. Finish Time 

This should match stop_time as specified in the define_cathode 
section of the code. 


c. Space Charge 

This option allows the user to turn the space charge on and off. If zero, 
the Poisson correction becomes a Laplace correction. The fields generated by the 
simulation will match those learned in an introductory course in electrostatics. 

With the space charge on, spiffe calculates the field from each 
macroparticle and its effect on the external field as well as neighboring particles. This 
requires considerably more computational steps and much longer run times. 

d. Sample Code 

The code used for both the box and parallel plate geometries is shown 
below. Finish_time always matched stop_time (from define_cathode 
section). Time is given in seconds. Section 4.19 of [14] provides a description of 
additional integration options and their default values. This section of code should be 
added at the bottom of the text file, after all the other code has been written. This is the 
only section of the code that is not optional. If this section is not included, spiff e will 
not calculate anything. Integration step size and run time were chosen using trial and 
error to find the best balance of accuracy and reasonable run times. 

Sintegrate 

dt_integration=le-14, 
start_time=0, 
finish_time=5e-10, 
space_charge=l, 


Send 
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G. 


POSTPROCESSING 


This section provides an introduction to the most commonly used post-processing 
tools from SDDS Toolkit. The SDDS Toolkit user’s guide [15] provides more detailed 
information and options. The intent is to guide the user to useful commands, which can 
then be learned in more detail from [15]. 

1. sddsquery 

This command is helpful for determining the type of data contained in any output 
file. Typing “sddsquery box.scO” at the command prompt will return a list of the data 
name, units, symbol, and type contained in the file box.scO. The data name is the 
terminology used to call on the data with other commands (e.g., sddsplot). 

Data is considered either a parameter or a column. The simplest way to explain 
this is through an example. The file box.scO was created from a define_screen 
command. This command captures beam data (r, p^, pr, t, q, etc.) at each time step. The 
output file is separated into tables; each table represents data from one time step. The 
tables are labeled with the time step of the data they contain; time step is the parameter. 
The beam data contained in the table is separated into columns. 

2. sdds2spreadsheet 

This command has been replaced 
understanding the layout of an output file, 
the following at the command prompt: 

sdds2spreadsheet box.scO 

The delimiter option (-delim=”l”) is 
data into other programs. 

3. sddscombine 

The output files are separated into tables, each with a set of parameters as the 


by sddsprintout, but it is helpful for 
To view box.scO as a text file (box.txt), type 

box.txt -delim="I" 

entirely optional. It is helpful for importing 


header. The user’s guide [15] describes using this command to combine several files into 
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a single file. Another option is to use it to combine the separate tables in one file into a 
single table. This is helpful because many programs will only see the first table during a 
data import. This command was the first step in calculating current through a screen. 
The following command combines the separate tables of box.scO into a single table in the 
file box.comb. 

sddscombine box.scO box.comb -merge -overwrite 

4. sddsprocess 

Spif f e data can be exported to a spreadsheet or Matlab® for processing, but the 
sddsprocess command performs many of the same functions. This command was 
used to calculate the current density through each screen. First, a screen was placed in 
the beam path using the def ine_screen command in the run file. Next, the data from 
the screen file (box.scO) was combined into the file box.comb. This file contains several 
columns of beam data, but the relevant columns for calculating current are time (t) and 
charge (q). Each time a particle passes through the screen, its charge and the time it 
passes through the screen are recorded. Calculating the current simply requires 
knowledge of the total charge that passed through the screen in a given amount of time. 
The current density is calculated by dividing by the area of the beam. This is all 
accomplished with one sddsprocess command. The output is sent to the file box.pro. 

sddsprocess box.comb box.pro -process=q,sum,qsum 
-process=t,spread,tspread -define=parameter,I,"qsum tspread 
/",units=A -define=parameter,J,"I 5e-3 sqr pi * 
/",units=A/m$a2. 

The notation “m$a2” equates to m^. The notation “$a” is telling the program to 
place the next character in the superscript position. The subscript position is 
accomplished using “$b”. The text is returned to the normal position using “$n”. 

Additional parameters may be defined by adding more -process options to the 
line above. To calculate the average longitudinal momentum (pz), add “- 
process=pz, average, pravg" to the line above. Average is a built-in process 
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for sddsprocess, pz is the name of the column data that is being processed, and 
pravg is the name of the new parameter as defined by the user. 

5. sddscollapse 

The current density (7) that was calculated using sddsprocess was defined as 
a parameter. If several parameters are defined (e.g, rmax, pxmin, pzmax, pzavg), it is 
helpful to use the sddscollapse command. This command will collapse all the 
parameters from the original file (e.g., box.pro) into a set of columns. The column data 
from the original file will not be transferred. The columns of the new file (box. coll) will 
be created from the parameters of the original file. 

sddscollapse box.pro box.coll 

6. sddsprintout 

The data contained in any sdds file can be viewed using sddsprintout. The 
user must specify which data to print; this can be columns, parameters, or both. To view 
the current density (J) and average longitudinal momentum (pzavg) from box.coll, type 
the following command at the command prompt. 

sddsprintout box.coll box.txt -col=(J,pzavg) 

The data is printed to the file box.txt, which the user must open to view the 
required data. 

7. sddscontour 

This command was used to create the electric field contour plots 
throughout this report. This is the simplest plotting tool available in SDDS Toolkit. A 
basic contour plot of can be obtained using the following command: 

sddcontour -quantity=Ez box.fields 
The number of contour lines may be adjusted by adding the “-contours” switch: 

sddscontour -quantity=Ez -contours=50 box.fields 


94 



The contour plot may be shaded by adding the “-shade” switch: 

sddscontour -quantity=Ez -shade box.fields 

8. sddsplot 

This command has a multitude of options available for the user. It can create 
quick plots for viewing on-screen or it can create publication-quality graphics. One of 
the most useful ways to utilize this command is to view beam snapshots. First, snapshots 
must be added to the run program using def ine_snapshot s. The plotting command 
must call on the particle location (z, r) from the snapshot file (box.snap) as well as the 
coordinates in the boundary file {box.dbnd). To obtain an on-screen plot, type the 
following command at the prompt: 

sddsplot -col=z,r -graph=dot -spllt=page -sep box.snap 
-col=z,r -graph=llne -spllt=page -sep -omnl box.dbnd 

Snapshots are viewed by pressing “n” for next and “p” for previous. If -sep is 
not used, the snapshots will be plotted on top of each other rather than in separate 
windows. This is helpful for showing the beam envelope for the entire simulation. 

A publication-quality graphic is obtained by requesting a white background and 
sending the plot to an output file {snap.png). 

sddsplot -col=z,r -graph=dot -spllt=page -sep box.snap 
-col=z,r -graph=llne -spllt=page -sep -omnl box.dbnd 

This command will only plot the first snapshot (first table of data). To plot the 
remaining snapshots, the snapshot file (box.snap) must be split into separate files using 
the sddsspllt command. 

9. sddssplit 

Output files contain several tables of data, with each table corresponding to a 
different time step. The sddsspllt command creates several new files from the 
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original file. Each new file contains data from one table of the original file. To split 
box.snap into several individual snapshot files, type the following command at the 
prompt: 

sddssplit box.snap -rootname=box -digits=3 - 
extension=snap 

If the original file contains 11 tables of data, this command will create 11 files 
named boxOOl.snap, box002.snap, ..., boxOll.snap. Each of these snapshot files may 
then be plotted to an output file using the sddsplot command described previously. 

H. TROUBLESHOOTING 

1. General Troubleshooting 

The boundary output should always be compared before proceeding with any 
simulations or troubleshooting. This was described in Section E.2.c of this Appendix. 

One of the easiest methods for verifying a working code is by viewing beam 
snapshots. If there are no particles in any of the snapshots, something is probably wrong 
with the code. The most likely causes are: 

• Voltage is not high enough to accelerate the particles in the given amount 
of simulation time. This is fixed by increasing either the voltage or 
simulation time. 

• Cathode stop_time does not match simulation finish_time. This 
only seemed to be a problem on Windows-based machines. 

• Cathode location might be inside a boundary. Make sure that the discrete 
boundary matches the ideal boundary and verify the location of the 
cathode. 

• Wrong dimensions. Verify all dimensions are in meters, including cavity 
boundaries and cathode location. 


96 



A second method of checking for particles is by viewing screen data. Use 
sddsprintout or sdds2spreadsheet to view the data. If the file is empty, no 
particles have passed through the screen in the direction indicated. It may be helpful to 
add another screen at the same location to record backward particles. 

Viewing the electric field contours is another way to verify that the code is 
working as expected. Run a simulation without particles to verify that the electric field is 
established as expected. Running the simulation with the space charge off is another 
option to check the electric fields. If the electric fields are not established as expected, 
verify the potentials and dimensions of the geometry file. 

2. Files not Found 

Reference [14] shows the use of both commas and semi-colons following a 
command/namelist. Use caution when using semi-colons. For example, using 


boundary_output=box.bnd, 
creates an output file named box.geo, but using 

boundary_output=box.geo; 

creates an output file name box.geo; (note that the semi-colon is now a part of the file 
name). 
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APPENDIX B SCREEN DATA 


Jreq (kA/lTI^) 

Screen 

Jdei (kA/m^) 

Pz,avg (rHeC) 

Pr.avg (rHeC) 

rmax(mm) 

5 

0 

4.991512 

6.259E-03 

0 

4.999995 


1 

4.990803 

1.399E-02 

0 

4.999995 


2 

4.991676 

1.877E-02 

0 

4.999995 

50 

0 

49.91512 

6.259E-03 

0 

4.999995 


1 

49.90803 

1.399E-03 

0 

4.999995 


2 

49.91676 

1.877E-02 

0 

4.999995 

100 

0 

99.83024 

6.259E-03 

0 

4.999995 


1 

99.81605 

1.399E-03 

0 

4.999995 


2 

99.83352 

1.877E-02 

0 

4.999995 

150 

0 

149.7454 

6.259E-03 

0 

4.999995 


1 

149.7241 

1.399E-02 

0 

4.999995 


2 

149.7503 

1.877E-03 

0 

4.999995 

200 

0 

199.6605 

6.259E-03 

0 

4.999995 


1 

199.6321 

1.399E-02 

0 

4.999995 


2 

199.667 

1.877E-02 

0 

4.999995 

250 

0 

249.5756 

6.259E-03 

0 

4.999995 


1 

249.5401 

1.399E-02 

0 

4.999995 


2 

249.5838 

1.877E-02 

0 

4.999995 

300 

0 

299.4907 

6.259E-03 

0 

4.999995 


1 

299.4482 

1.399E-02 

0 

4.999995 


2 

299.5006 

1.877E-02 

0 

4.999995 

350 

0 

349.4058 

6.259E-03 

0 

4.999995 


1 

349.3562 

1.399E-02 

0 

4.999995 


2 

349.4173 

1.877E-02 

0 

4.999995 

400 

0 

399.3209 

6.259E-03 

0 

4.999995 


1 

399.2642 

1.399E-02 

0 

4.999995 


2 

399.3341 

1.877E-02 

0 

4.999995 

450 

0 

449.2361 

6.259E-03 

0 

4.999995 


1 

449.1722 

1.399E-02 

0 

4.999995 


2 

449.2508 

1.877E-02 

0 

4.999995 

500 

0 

499.1512 

6.259E-03 

0 

4.999995 


1 

499.0803 

1.399E-02 

0 

4.999995 


2 

499.1676 

1.877E-02 

0 

4.999995 


Table 27. Screen data, box, space charge off. 
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Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz.avg (rngC) 

Pr,avg (rHeC) 

rmax(mm) 

5 

0 

4.990588 

6.229E-03 

-1.658E-08 

5.000001 


1 

4.989335 

1.396E-02 

-9.347E-09 

5.000076 


2 

4.990293 

1.876E-02 

-1.141E-07 

5.000168 

50 

0 

49.83783 

5.976E-03 

-1.442E-07 

5.000058 


1 

49.78213 

1.372E-02 

1.204E-09 

5.000835 


2 

49.78434 

1.871E-02 

1.028E-06 

5.001752 

100 

0 

99.51599 

5.661E-03 

-2.198E-07 

5.000131 


1 

99.31531 

1.344E-02 

3.697E-07 

5.001751 


2 

99.28589 

1.864E-02 

-1.549E-06 

5.003623 

150 

0 

148.9618 

5.293E-03 

-1.830E-07 

5.000217 


1 

148.4779 

1.313E-02 

1.378E-06 

5.002786 


2 

148.3694 

1.857E-02 

-1.314E-06 

5.005690 

200 

0 

197.9707 

4.840E-03 

1.166E-07 

5.000320 


1 

196.9827 

1.278E-02 

3.378E-06 

5.004042 


2 

196.7009 

1.850E-02 

1.782E-07 

5.008131 

250 

0 

245.6037 

4.211E-03 

9.568E-07 

5.000465 


1 

243.7504 

1.239E-02 

7.391E-06 

5.005637 


2 

243.2093 

1.842E-02 

4.033E-06 

5.011119 

300 

0 

270.216 

3.484E-03 

3.942E-06 

5.000716 


1 

267.0417 

1.211E-02 

1.445E-05 

5.007847 


2 

266.0374 

1.838E-02 

1.161E-05 

5.015166 

350 

0 

271.2137 

3.385E-03 

3.364E-06 

5.001017 


1 

268.0365 

1.209E-02 

1.595E-05 

5.010126 


2 

266.778 

1.838E-02 

1.373E-05 

5.018893 

400 

0 

270.7643 

3.343E-03 

4.043E-06 

5.001373 


1 

267.4595 

1.208E-02 

1.700E-05 

5.011907 


2 

266.4258 

1.838E-02 

1.453E-05 

5.021549 

450 

0 

271.6229 

3.326E-03 

5.416E-06 

5.001679 


1 

267.3812 

1.208E-02 

1.837E-05 

5.012804 


2 

267.0899 

1.838E-02 

1.565E-05 

5.022751 

500 

0 

271.5669 

3.314E-03 

4.669E-06 

5.002326 


1 

268.3792 

1.208E-02 

1.869E-05 

5.015046 


2 

266.9832 

1.838E-02 

1.628E-05 

5.025171 


Table 28. Screen data, box, space charge on. 
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Jreq (kA/lTI^) 

Screen 

Jdei (kA/m^) 

Pz,avg (rHeC) 

Pr,avg (rHeC) 

rmax(mm) 

5 

0 

4.991512 

6.256E-03 

0 

4.999995 


1 

4.990803 

1.399E-02 

0 

4.999995 


2 

4.991676 

1.877E-02 

0 

4.999995 

50 

0 

49.91512 

6.256E-03 

0 

4.999995 


1 

49.90803 

1.399E-02 

0 

4.999995 


2 

49.91676 

1.877E-02 

0 

4.999995 

100 

0 

99.83024 

6.256E-03 

0 

4.999995 


1 

99.81605 

1.399E-02 

0 

4.999995 


2 

99.83352 

1.877E-02 

0 

4.999995 

150 

0 

149.7454 

6.256E-03 

0 

4.999995 


1 

149.7241 

1.399E-02 

0 

4.999995 


2 

149.7503 

1.877E-02 

0 

4.999995 

200 

0 

199.6605 

6.256E-03 

0 

4.999995 


1 

199.6321 

1.399E-02 

0 

4.999995 


2 

199.667 

1.877E-02 

0 

4.999995 

250 

0 

249.5756 

6.256E-03 

0 

4.999995 


1 

249.5401 

1.399E-02 

0 

4.999995 


2 

249.5838 

1.877E-02 

0 

4.999995 

300 

0 

299.4907 

6.256E-03 

0 

4.999995 


1 

299.4482 

1.399E-02 

0 

4.999995 


2 

299.5006 

1.877E-02 

0 

4.999995 

350 

0 

349.4058 

6.256E-03 

0 

4.999995 


1 

349.3562 

1.399E-02 

0 

4.999995 


2 

349.4173 

1.877E-02 

0 

4.999995 

400 

0 

399.3209 

6.256E-03 

0 

4.999995 


1 

399.2642 

1.399E-02 

0 

4.999995 


2 

399.3341 

1.877E-02 

0 

4.999995 

450 

0 

449.2361 

6.256E-03 

0 

4.999995 


1 

449.1722 

1.399E-02 

0 

4.999995 


2 

449.2508 

1.877E-02 

0 

4.999995 

500 

0 

499.1512 

6.256E-03 

0 

4.999995 


1 

499.0803 

1.399E-02 

0 

4.999995 


2 

499.1676 

1.877E-02 

0 

4.999995 


Table 29. Screen data, parallel plate, space charge off. 
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Jreq (kA/m^) 

Screen 

Jdei (kA/m^) 

Pz.avg (rHeC) 

Pr,avg (rHeC) 

rmax(mm) 

5 

0 

4.990464 

6.213E-03 

-1.186E-08 

5.000000 


1 

4.988699 

1.393E-02 

2.341E-10 

5.000071 


2 

4.989187 

1.871E-02 

-9.912E-08 

5.000159 

50 

0 

49.7854 

5.793E-03 

-6.424E-08 

5.000056 


1 

49.64779 

1.333E-02 

3.014E-07 

5.000874 


2 

49.57408 

1.819E-02 

-5.570E-07 

5.001852 

100 

0 

99.08017 

5.239E-03 

6.545E-08 

5.000176 


1 

98.39032 

1.261E-02 

1.659E-06 

5.002332 


2 

98.02934 

1.760E-02 

3.218E-07 

5.004700 

150 

0 

146.5805 

4.517E-03 

6.389E-07 

5.000421 


1 

144.4719 

1.182E-02 

5.052E-06 

5.005229 


2 

143.7236 

1.702E-02 

3.715E-06 

5.009518 

200 

0 

176.3814 

3.712E-03 

2.529E-06 

5.001135 


1 

174.034 

1.118E-02 

1.137E-05 

5.009485 


2 

174.101 

1.660E-02 

1.071E-05 

5.017878 

250 

0 

184.8474 

3.284E-03 

4.302E-06 

5.001837 


1 

183.1321 

1.091E-02 

1.608E-05 

5.013340 


2 

183.5346 

1.645E-02 

1.582E-05 

5.023391 

300 

0 

187.3599 

3.112E-03 

4.479E-06 

5.002438 


1 

184.9992 

1.083E-02 

1.831E-05 

5.015919 


2 

185.8662 

1.641E-02 

1.830E-05 

5.027015 

350 

0 

188.0085 

3.050E-03 

4.414E-06 

5.003105 


1 

185.6726 

1.080E-02 

1.957E-05 

5.017951 


2 

186.3659 

1.640E-02 

1.974E-05 

5.029710 

400 

0 

187.8295 

3.028E-03 

9.802E-06 

5.003967 


1 

185.6414 

1.080E-02 

2.188E-05 

5.020319 


2 

186.3309 

1.640E-02 

2.123E-05 

5.032293 

450 

0 

187.9653 

3.013E-03 

8.929E-06 

5.004827 


1 

185.9714 

1.079E-02 

2.201E-05 

5.020169 


2 

186.7722 

1.639E-02 

2.185E-05 

5.032424 

500 

0 

187.8994 

2.996E-03 

6.467E-06 

5.006683 


1 

186.4457 

1.077E-02 

2.156E-05 

5.024437 


2 

187.3612 

1.638E-02 

2.210E-05 

5.037074 


Table 30. Screen data, parallel plate, space charge on. 
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